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Introduction 


When we were solving sets of simultaneous equations like 
2x, + 3x, — 4x, =5 
X,+2x, + x; =3 
3x, — x, + 4x; =6, 


we found we were particularly interested in the rectangular blocks of numbers 


2 3°-4 5 
1 2 1| and |3 
A Se 6 


These rectangular blocks are called matrices (plural of matrix). There are other 
situations where such blocks of numbers occur. For instance, look at the 
information in Table 1 below showing the prices of different articles in competing 
shops. (Prices are given in pence.) 


Table 1 
Dessert Dessert Small 
Fork spoon fork knife 
52 48 48 58 
pa! 45 j 45 if 37 
| 54 47 ar 49 | 56 


60 52 48 48 58) 


63 54 47 49 56] 
You will meet a matrix like this in Section 2, 
Another example of a simple matrix can be obtained from the vector 
; t= di + aaj + ayk. 


The essential information about the magnitude and direction of r is contained in 
the three numbers a,, a and ay. So it is quite usual to represent r by a column 
ay 
matrix a, |, or sometimes by a row matrix [a, a, a3). 
as 


It is because vectors are often represented by such column (or row) matrices, that 
such matrices are in fact usually referred to as column yectors (or row vectors). 
‘These terms are extended to refer to any matrix consisting of a single column (or 
row) of any length. For example, [1 24 4 7] is referred to as a row 

vector. 


Manipulation of matrices like the ones above occurs sufficiently often for it to be 
worth while to develop an algebra of matrices which will enable us to refer simply 
and briefly to the manipulations we carry out, 


You saw in the unit on complex numbers how it was possible to create an algebra 
for things other than real numbers, with suitable definitions for *+° and *x*. We 
can manipulate complex numbers in the same way as real numbers, simply 
because they obey the same laws. When we have defined what we mean by matrix 
‘addition’ and ‘multiplication’, we will find that the rules for matrix algebra are 
largely similar to those for real and complex numbers. The main purpose of this 
unit is to investigate these and to learn how to use the resulting algebra. For 
comparison, the main laws for real numbers are given in the Handbook, 
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Study guide 


You must make sure that you understand the contents of Sections 1 and 2 and 
Subsection 3.1 (the tape subsection) before you watch the television programme. 
This is probably about three hours’ work. 


If you have studied matrices before, you should find this a reasonably easy unit. It 
is for this reason that there are not too many matrix multiplication examples 
included in the body of the text. However, if you have never met matrix 
multiplication before, you will probably find that you need to do some more 
examples, as it takes a bit of time to get facility with this operation. It is for this 
reason that there are plenty of multiplication examples at the end of Section 2. In 
general, the end of section exercises should be used, either as extra practice if you 
need it, or for revision purposes. The end of unit test is there for you to see if you 
have assimilated the contents of the unit! 


It is expected that you will study Sections | and 2 in the first working session. 
This will be a bit heavy if you haven't studied matrices before, but some of the 
later sections will take much less time. 


Finally, Section 5 is about determinants. Determinants are not matrices, but to 
study this section you need to know the definition of a matrix. 


1 Some simple matrix operations 


1.1 Matrix notation 


You will remember from Unit 9 on simultaneous equations that a matrix is a 
rectangular block of numbers. When we wanted to refer to a general set of two 
equations in two unknowns, we wrote 

aX, + b,x, =c, 

a,x, + b,x, =¢3, 


and when we wished to use the matrix of the left-hand side coefficients we wrote 


This notation is fine when discussing small sets of equations or small matrices, but 
we would run out of letters of the alphabet if we wished to talk about larger 
matrices. It is for this reason that we adopt the notation 

ausX, + 442%, =, 

Gy4X1 + Gz2X, = by. 
This very convenient notation means that we never run out of names for the 


coefficients, and can cope with any number of equations or any size matrix in this 
way. 


It is usual practice to use a capital letter when referring 10 matrices. In this course, 
we shall use a capital letter in bold type. Thus the matrix 


bes al 
ay ay, 


will be referred to briefly as A. Other matrices could be called X or B, for example. 
The only exception to this notation is that used for matrices consisting of a single 
column. It is common practice to refer to these ‘column vectors’ using a lower-case 
Jetter. We shall use a lower-case letter in bold type for these. Thus the column 


vector 
*y 
| 


will be referred to briefly as x. 
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Now, when we use the letter A to represent a matrix, it tells us nothing about its 
size or its contents. Suppose we wish to refer to a matrix which has two rows and 
three columns. We call this a two by three matrix and write it as ‘2 x 3’, 


Note: it is important to specify the number of rows before the number of columns. 


Exercise 1 

Write down examples of matrices of the following size, 
(i) 2x2 (i) 3x4 (iii) 4x1 (iv) 1x 2. 
[Solution on p. 55} 


Each individual entry in a matrix is called an element, The element qy in a matrix 
Q refers to the element in the ith row and jth column. 


Example 1 
The elements of a 2 x 4 matrix Q are labelled in the following way: 


q-(% M2 ds td 
424 qaa 425 q24 


In this way, each individual element in a given matrix has its own name, and we 
know where to find it. 


Exercise 2 
How would we refer to the element in the fourth row and fifth column of a matrix P? 
[Solution on p. 55} 


In general, a matrix A with m rows and n columns is called an 
m x n matrix: 


Sie Sie pee ain 

a a gy, 
A=|% 22 Falk 

Am Oma Ms mn 


In this matrix, a,, is the element in the ith row and jth column. 


In the case of a column vector x, it is only necessary to specify one subscript. For 
instance, we know that x, is the element in the third row of a column vector x. 


1.2. Equality of matrices 


We define ‘A = B’ to mean that A and B are identical, ie. 
(a) that the matrices A and B are the same size, 
(b) that the corresponding elements are equal. 


So, if am |e mia a and B 


a4 S00 as 


then A = B means that 


4,=1, ay=2, 


a, = 0, a= 4, 


In this instance, the single matrix equation provides us with a shorthand notation 
for writing down six real-number equations. 
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Exercise 3 
Write down the real-number equations implied by the statement A = B, if: 


i) re é| and s-| § a 


ax + by c 
aas(e2?) w 2~(4] 


[Solution on p. 55) 


13 Matrix addition 


You know from Unit 14 on yector algebra, that the sum of the two vectors 
a=aitajt+ak and b=bi+b,j+b sk 


a+b=(a, +b, ji + (a, + by)j + (ay + b5)k. 
In our matrix notation, this would be written as 


ay by a, +b, 
a,| +] b,| =| a, +6, }. 
ay bs a, + bs 


So, if this matrix representation of vectors is to be at all useful, the definition of 
matrix addition would have to include this statement above. A fairly intuitive 
extension of this gives us the following definition. 


If A and B are two matrices of the same size, we define 
A+B 


to mean the single matrix formed by adding the corresponding 
elements of A and B. 


Example 2 
im 3 1 0 241 3+0 
1 4/+]/0 4/=] 140 444 
- 5 1 “) ed 5+3 
3 3 
= |1 8 
0 8 


A= (1 2) and 


we would not be able to add these matrices as, by definition, matrices must be the 
same size before we can add them. 


Another example of matrix addition was seen in the Gaussian elimination method 
in Unit 9. Each row of the matrix we used in the elimination process can itself be 
considered as a matrix. We were quite cheerfully adding one row to another so as 
to obtain a new row during each stage of the elimination process. 


Exercise 4 
Given that 
2 ‘| =f 4 0 Fy 3 
A=|4 3 B= 1 1 0 C=/1 
1 5 0 1 1 6 3 


3 

2|, 

6 
aye cles a 
sd VG a VS eens [he 
find (where possible): 
(i) A+B, (ii) A+ C, (iii) C+ x, (iv) (x+y) +z 
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Exercise 5 


Write down the real-number equations implied by x + y =z, where x, y and z are the 
matrices defined in Exercise 4. 


Exercise 6 

A matrix C is formed by adding two m x n matrices A and B. 
(i) What size will C be? 

(ii) How do we form the element Cy? 

[Solutions to Exercises 4-6 on p. 55] 


Before we can do any algebra with our matrices, we must check whether the rules 
for addition which operate for real numbers are true for matrices. In the following 
proofs we shall assume that the matrices are of the same size. 


1. To prove the commutative law for matrices for addition: 
A+B=B+A. 


Proof 


A typical element in the ith row and jth column of A + B is a, + by. The 
equivalent element for B + A is b,, + a,,. 


So, as we know from the commutative law for real numbers that 
4, + by = bj; + ay, it follows that 


A+B=B+A. 
2. To prove the associative law for matrices for addition: 
A+(B+C)=(A+B)+C. 


Proof 


A ‘typical’ element of A + (B + C) is a,, + (b,, + ¢,)). The equivalent element in 
(A + B) + C is (a,, + b,) + ¢,. Now, from the associative law for real 
numbers, we know that 


yy + (by + Cy) = (ayy + by) + cy. 
Hence we derive the associative rule for matrices: 
A+(B+C)=(A+B)+C. 
Having proved these two laws for matrices, we can now use them when required. 


It is useful at this stage to define the zero matrix 0, and the negative of a matrix A. 
These are analogous to the number zero, and the negative of a number a, in our 
ordinary number system. 


We define a zero matrix as a matrix consisting entirely of zeros. 
Such a matrix (whatever its size) is denoted by 0, and has the 
property that 


A+0=A. 


The size of a zero matrix is fixed by the context in which it arises. For instance, if 
A above were an m x n matrix, then the ‘0° would also be an m x n matrix, 


A matrix —A is one whose elements are the negatives of those of 
A. It has the property that 


A+(-A)=0. 


For instance, if A = [1 —2], then —A = [-1 2). 


Whenever we write A — B, this will be taken to mean A + (—B). In practice we 
just subtract the corresponding elements of B from A. 
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Exercise 7 


Given that A ; a i} write down the matrix —A, and the zero matrix such that 
A+(-A)=0. 
Exercise 8 
ay 2 ( 
WA=|1 2/andB=|-1 1|,findaA-B. 
Sound 2 - 


[Solutions to Exercises 7 and 8 on p. 55) 


1.4 Multiplication of a matrix by a number 
The effect of performing the operation 
A+A 


fs to obtain a matrix whose elements are double those of A. Similarly, if we add 
A +A + A, we obtain a matrix whose elements are three times those of A. It is 
tempting to think of A + A as 2A, and A + A + A as 3A, and we extend this idea 
to make the following definition. 


If k is a real number, and A is any matrix, then we define 
kA (or Ak) 


to be a matrix (which is the same size as A) in which each 
element is the corresponding element of A multiplied by k. 


It should be noted from this definition that 
(-I)A=-A and 0A =0. 


Example 3 
2 1 3 10 3s 15 
wa-[j 2 Hi then sa=[75 10 Al 


With this definition of multiplication of a matrix by a number, it can be shown 
that: 


(i) (ky ka)A =k, (kA), 

(ii) k(A + B)=kA + kB, 

(iii) (ky +k )A=kA+ kA. 

These laws, together with the laws for matrix addition, mean that, provided our 
matrices are of the same size and we confine ourselves to the algebra in this 


section, we can manipulate matrix algebra in much the same way as real number 
algebra. 


Exercise 9 


2 1 sil 0 

Given k =5,m=2,A=|3 2) and B= 1 | write as a single matrix the 
4 th 

following expressions: 


fi) kA, 
(ii) k(A + mB), 
(ii) kA — mB, 


(iv) mA + mB + kA — kB, 
[Solution on p. 55} 


Example 4 
We can find the matrix X from the matrix equation 
2(A+3X)-A=B+A 
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by using the laws and rules of this section. The references are to the list of rules 
given in the summary of this section. 


2A+6X-A=B+A (Summary 3(iv) then (iii) 
6X+2A—-A=B+A (Summary 3(i)) 
6X+A =B+A. (Summary 3(v)) 
Adding —A to both sides of the equation, we get 
6X=B (Summary 3(v)) 
and dividing both sides by 6, we get 
X=1B, (Summary 3(iii)) 


In fact, as you can see, this solution involves only the operations discussed in this 
section. We work with this equation in exactly the same way as we would with an 
equation involving numbers. 


Exercise 10 
(i) Simplify the expression 
2(A + 3B) — 4(A — B). 
(ii) Solve the following equation for X: 
X+3(A + X)=2X—A. 
[Solution on p. 55] 


Summary of Section 1 


1, Notation 


(i) A matrix is represented by a capital letter in bold type, eg AL BS. 
(ii) Column vectors are represented by lower case letters in bold type, e.g. a,b... 
(ili) A matrix.A with m rows and n columns is called an m x n matrix, 

(iv) aj, is the element in the ith row and jth column of a matrix A. 


2. Algebra 


(i) Equality: A = B means that A and B are the same size, and their 
corresponding elements are equal. 


(ii) Matrix addition: if two matrices A and B are the same size, then 
A+B 
is the single matrix formed by adding the corresponding elements of A and B. 
(iii) kA is the matrix whose elements are k times those of A. 
Two particular matrices are defined: 
(iv) a matrix 0 (consisting entirely of zeros) with the property that A +0 = A; 


(v) a matrix —A whose elements are the negatives of those of A. It has the 
property that A + (—A)=0. 


3. Laws 
Given matrices A, B, C (all the same size), then: 
(i) A+B=B+A, 


(ii) A+(B+C)=(A+B)+G 
(iii) (kyk.)A = k,(k2A), 

(iv) k(A+B)=kA + kB, 

(v) (ky +k,)A =k,A +k A. 


In this way, we start to build up an algebra for matrices very similar to the 
algebra for real numbers. 
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End of section exercise 


Exercise 11 

2 4 2 x 2 
peal al i} x-[eh oh 
(i) Find, if possible, 


(a) A+B, (b) B-C, (c) C + 2x (d) x + 2y, (e) A— 2B+C. 


(ii) Find the matrix P if 
2(P + 2A) + B= 3C. 
[Solution on p. 55] 


2 Matrix multiplication 


2.1 Motivation 


Example 1 

milk bread sugar 

(pints) (loaves) (kg) A B 
John 6 3 2 milk 18 17 
Jane 10 2 2 bread 30 32 
Joyce “4 | 2 1 sugar 35 33 
Jim 15 3 Table 2 

Table 1 


Table 1 above shows four people's weekly milk, bread and sugar requirements, and 
Table 2 shows the cost (in pence) of these items in two supermarkets A and B. 
There is another table which can be constructed from these two tables—a table 
which shows how much the four people will have to pay for their requirements in 
either supermarket. 


Table 3 


To work out how much it would cost John to shop in A, he buys 
6 pints of milk at 18p each, 
3 loaves at 30p each, and 
2kg sugar at 35p each, 
giving a total cost of 
(6 x 18) + (3 x 30) + (2 x 35) 
= 268p. 
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To work out how much it would cost Joyce to shop in B, we take Joyce's row in 
Table 1 (her requirements), and the ‘B’ column in Table 2 (the cost of these 
requirements), and obtain the cost 

(7 x\17) +,(2 x 32) + (1 x 33) 

= 216p. 
Exercise 1 
Fill in the rest of the entries in Table 3, 
[Solution on p. 56} 


This combination of the two tables could be expressed in matrix terms: 


6 3 2 18 17 268 264 
10 2 2] combined with 130 32) gives | 310 300 
Te iat el 3533 221-216] 
is) 495 482 


This process, of combining rows of the left-hand matrix with columns of the right- 
hand matrix to produce elements of a new matrix, is the operation which we call 
matrix multiplication, This is probably the operation which we perform most 
frequently when using matrices, You will see some examples of its use in this unit, 


t 


a,b, + aybz + asby 


2.2 The definition of matrix multiplication 
When we write 
C= AB, 


we mean we are combining (multiplying) two matrices A and B to obtain a new 
matrix C. Each element of this matrix is formed by taking a row of the left-hand 
matrix A, and combining it with a column of the right-hand matrix B, along the 
lines indicated in the previous section. The following example should show how 
the individual elements of such a matrix C are formed. 


Example 2 
Find AB if 
6 =! 
A= [; : 4 and B= 0 2 
=I 1 
Solution 
Let 


wwac-[ 3] 


Cay 22 

where ¢,, is formed by taking row 1 of A and combining it with column | of B, 
giving 

ey, =(2 x 6) + (1 x0) + (3 x -1)=9, 
To get c,,, combine row 1 of A and column 2 of B, to give 

ey. =(2x -1) 4+ (1 x 2)4+(3x 1) =3, 
To get c,,, combine row 2 of A and column 1 of B, to give 

Cn, = (4 x 6) + (5 x 0) + (6 x —1)=18, 
To get c,,, combine row 2 of A and column 2 of B, to give 

32 = (4x -1)+ (5x 2)4+(6 x 1)=12, 


You will notice the analogy 
with the formula for 
calculating scalar products of 
vectors. 
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Thus we obtain 


| ea 
Pare) 
ee 
aw 
ey 
—0 
- ODN 

i 
ples 
eed 
W 
=— 
cay 
Bw 
cee 


A B = c 


So, to find an element ¢,, of a matrix C = AB, take row i of the left-hand matrix A, 
and combine it with column j of the right-hand matrix B in the manner described 
above. 

Exercise 2 

Find AB if 


[Solution on p. 56) 


There are more exercises on matrix multiplication later in this section. Make 
certain that you do enough of them to feel confident that you can do the process. 


The method above actually determines the size of C = AB. For instance in 
Example 2, to obtain an element c;,, we would have had to combine row 3 of A 
with row 1 of B. As A has no third row, neither has C, So we conclude that C has 
the same number of rows as the left-hand matrix A. Similarly, to obtain an element 
¢, 3, we would have to combine row 1 of A with column 3 of B, As B has no third 
column, neither has C. So we conclude that C has the same number of columns as 
the right-hand matrix B. 


Example 3 


If A is a 2 x 3 matrix, and B is a 3 x 4 matrix, then AB will be a 2 x 4 matrix. 
For instance, 


Exercise 3 
Find the size of AB, if A is a 3 x 5 matrix, and B is a 5 x 2 matrix. 
[Solution on p. 56) 


Finally, the rule for combining the elements of A and B has automatically imposed 
a restriction on the size of matrices we can multiply together. For instance, if a 

b, 
typical row of A were [a, a, a3; 4,4] and a typical column of B were | b, |, 

b; 
then we wouldn't be able to do the specified row-column operation, as we would 
have an element a, left over. So we have the natural restriction, that to be able to 
Jorm AB, the number of columns of A must equal the number of rows of B. 


Example 4 


If A is a 2 x 3 matrix, then if AB is to be formed, B will have to have 3 rows. For 
instance, the products 
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can all be formed, whereas to form the product 


| 


is not possible. 


(Note, there is no restriction on the number of columns that B may have, or the 
number of rows that A may have.) 


Exercise 4 
If A is a 2 x 4 matrix, then what size must B be to be able to form 
{i) AB (ii) BA? 


[Solution on p. 56} 


The restrictions on the sizes can be summed up as follows. 


Exercise 5 
Given the matrices 


(i) state which matrix products can be formed, 
(ii) state the size of the resulting matrix products. 


Exercise 6 
Form all the possible matrix products from the matrices 


ap ed ae of 


Exercise 7 
Given the matrices 
2 0 1 1 2-1 
A= |-1 2 1 and B=/-1 2 ABs 
3 4 3 Yae aa | 


find the matrix products AB and BA. 
[Solution to Exercises 5-7 on p. 56] 


The ideas of this subsection can be summarized in the following definition of 
matrix multiplication. 


The matrix C = AB can be formed if the number of columns of 
the left-hand matrix A equals the number of rows of the right- 
hand matrix B. The matrix C will have the same number of rows 
as A, and the same number of columns as B, and its elements are 
given by 


€4j = Gixb1y + Giada) + °° + indy, 
n 
= 2 dubs 
At 


where n is the number of columns of A. 
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2.3. The algebraic laws governing matrix multiplication 


Before we can freely use matrix multiplication, we need to know whether matrix 
multiplication is commutative and associative, and whether there is a distributive 
law. 


As you can see from the results of Exercise 7, AB and BA are not the same. Any 
counter-example would have been adequate to show that, in general, matrix 
multiplication is not commutative, that is, 


AB is not usually equal to BA. 


In fact, we could have seen this without doing an example at all. If we multiply an 
m x n matrix A to an n x p matrix B, we obtain an m x p matrix C. However, it is 
not even possible to form the product BA unless p is equal to m, 


As this is the first thing that we have discovered about matrices which differs 
significantly from real number algebra, we shall have to be very careful. 


For instance, if 
A=B (1) 


and we wish to multiply both sides of this matrix equation by another matrix C, 
then we have to decide either to multiply both sides of (1) by C on the left, to 
obtain 


CA = CB, (2) 
or we could (size permitting) multiply both sides of (1) on the right, to obtain 
AC = BC. (3) 


Both the matrix equations (2) and (3) would be logical consequences of (1), 
whereas there would be no justification in assuming that that AC and CB were 
equal. 

Despite this bad beginning to the investigation of the laws governing matrix 
multiplication, in fact the associative and distributive laws are exactly the same as 
those in real number algebra. However, as matrix multiplication is not 
commutative, there are two versions of the distributive law to prove: 


A(B + C) = AB + AC 

and 
(B+ C)A = BA +CA. 

The proofs of these, and of the associative law 
A(BC) = (AB)C 


are given as answers to the following exercise. They are there for completeness, but 
please don’t attempt to do them unless you are very brave, and have time on your 
hands! Exercises 16 and 17 at the end of the section will demonstrate these laws 
for three particular matrices. You are strongly advised to do these exercises now if 
you haven't had much practice at matrix multiplication before. 


Exercise 8 (Optional) 


Assume (for simplicity) that A, B and C are all n x n square matrices. Prove, using the 
definition of matrix multiplication given in the box on p. 14, that: 


(i) A(B+C) = AB+ AC, 

(ii) (B+ C)A =BA +CA, 

(iii) A(BC) = (AB)C. 

[Solution on p. 56) 

These laws, along with the laws of Section 1 and the inverse matrix in Section 4 
(the matrix equivalent of the number a~'), will enable us to do algebra with 


A consequence of the 
associative law, 

A(BC) = (AB)C, is that we 
can omit the brackets where 
appropriate and write ABC 
without fear of ambiguity. 
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matrices in much the same kind of way as we do with real numbers. The only 
difference will be that we shall have to be careful not to interchange expressions 
like AB and BA. 


2.4 The product of a matrix and a column vector 


The product of a matrix and a column vector produces another column vector— 
for example, 


2 7 1 1 15 
3 4 2 2] = | 9]. 
5 0 =1) [=1 6. 


Exercise 9 


2 1 2 
ra=[3 i ]nda=[7 sin Aq. 


[Solution on p. 57) 


Example 5 


3 1 x 
A= and x =|" "|, then 
4 2 x, 
3x, + x, 
Ax = ; lh 
4x, + 2x, 
In this last example, the column vector formed looks very like the left-hand side of 
a pair of simultaneous equations. It can be seen that any matrix A multiplied by 
an ‘unknown’ column vector x provides us with a column vector whose elements 


are linear combinations of the elements of x. This enables us to write any set of 
linear simultaneous equations in the form 


Ax =b. 
Example 6 


Given that A= [; | x= ie b= [*} 


write down the simultaneous equations represented by 
Ax =b. 
Solution 
1 -1][x,]_[o 
2 Sail vs 
gives us 
x,—x; ]_[0 
2xji+ 3x, J] 5] 


By definition of equality of matrices, this gives us the simultaneous equations 


x,- *,=0 
2x, + 3x, =5. 
Exercise 10 


Write down the simultaneous equations represented by Ax = b 


3) ft <2 x 2 
ffA=]1 2 1), x=|x,], b=] 1 
fe 0 ee -1 


[Solution on p. 57) 
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Example 7 
Express the equations 
2x, + x, =3 
x, —3x,=5 
in the form Ax = b. 
Solution 


Reversing the procedure of Example 6, the left-hand side of these equations can be 
written as 


2 1][ x, 
1 -3][x, J 
So the equations above can be written as 


2. Aqfs)|_[3 
1 -3]Lx,}]° [5 
which has the required form Ax = b. 


Exercise 11 

Write the equations 
2x,—3x,+ xy= 1 
Ixy + X,y-— X3= 2 
9x, + 3x, + 5x3=-1 

in the form Ax = b. 

[Solution on p, 57] 


The reason that this notation was not used in the unit on simultaneous equations 
was that the operation of matrix multiplication had not been discussed at that 
time. 

There is another context in which linear combinations and matrices arise. Suppose 
that we have two pairs of numbers (x,,y,) and (x,,2) such that 

Xp = X, +2y, 

} «) 
Y= 4x+ Vi 
These equations allow us to put values of x, and y, in the right-hand side of (1), 
and hence calculate the values of x, and y,. In a sense, we have ‘transformed’ 
(x,,y,) into (x2,)2), So a set of equations, like (1) above, is called a 
transformation—in this case a linear transformation. All linear transformations 
may be written in matrix form. For instance, the Equations (1) above can be 
written as 


El os 


Thinking of transformations in this matrix form enables us to manipulate them 
more easily. For instance, suppose we know that 


2 = si X3 = 3x, - * 
1 and (2), 
Y= 4x4 V1 2) Ya= X2+2y, e 


and we wish to express x, and y, in terms of x, and y,. 


We know that (1) and (2) above can be expressed in matrix form as 
x.]_f1 9 2]fxy x3/_ [3 -1][x, 
Pll lo = Cl-b al] « 


Then we can replace the expression for the column vector [=] from (3) into (4), 
2 


Note that A is the matrix of 
the left-hand side coefficients, 
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ee HE ES 
OE Ie 


using the associative law for matrices. Thus x, and y, can be expressed more 
simply in terms of x, and y, as 


BCs aE] 

3 9 4 

Rather more generally, suppose we perform a sequence of linear transformations of 
this kind: 


m=Ar, (3) ae 
=Br (4) bill 
Then by substituting the expression for r, in (3) into (4), we get 
ty; = B(Ar,) 
= (BA)r,, 


and hence we obtain expressions for x, and y, in terms of x, and y,. 

You will need this concept of combining successive matrix transformations using 
the multiplication process to understand what is happening in the television 
programme. 


Exercise 12 


1" 10 Sead =A x 
wal) i} s-[) | c-[ 1 Hl en ee -[r} 
express the following succession of matrix transformations: 
r,=Ar,, ,=Br, and r,=Cr, 
as a single matrix transformation 
14 =Qr,. 
Find the matrix Q, first in terms of A, B, C and then as a matrix with numerical elements, 
[Solution on p. 57] 


2.5 The matrix transpose 


The transpose of a matrix A (written as A’ in this course, or 
sometimes as A’) is simply the matrix whose columns are the 
rows of A. 


For instance, if 


You can see that the first row of A becomes the first column of A’ and the second 
tow of A becomes the second column of A’, 


Exercise 13 


2 3 3 
wal 0 2 


SS 
a 
J 
! 
! 
eon 
SME 


write down AT and B7. 
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Exercise 14 
Show for the matrices in Exercise 13 that 
(AB)? = BTAT. 
[Solutions to Exercises 13 and 14 on p. 57] 
The result 
(AB)" = BTAT 3) 


is always true, although it is not proved in this unit. We shall assume this result 
and show that a similar result is true for three matrices, 


Example 8 

To show that (ABC)™ = C™BTA’. 

Proof (ABC)? = (A(BC))” 
= (BC)"AT using (5) 
= C'BTAT using (5) 


This is just another example of how to use matrix algebra. 
The general result is 


(AyA5.:0A,)” ="ATAT |... ATA? 


An example of the use of a matrix transpose is in the theoretical matrix 
description of the linear programming problem in Unit 10. A typical linear 
programming problem would be: 

Minimize c,x, + ¢.x2 + ¢3X 5 

subject to the conditions 


yy X, + Gy 2X2 sb, 
4X, + az5X, Sb, 
X1,X2,X3 20. 


In matrix terms, this can be expressed as: 
Minimize [c, c¢, ¢3][x, 
*2 
X3. 


subject to the conditions 


i 4. 0 | * by 
a, 0 423] |x,| < | ba}, 
x3 by 


Xy,X2)X3 20. 
Rather more briefly, we say: 
Minimize c7x 
subject to the conditions 
Ax <b 
and x20. 
We are forced to use the notation ¢7 as, to save confusion, we stated at the 


beginning of the unit that the notation ¢ refers only to a column vector. So if we 
need a row vector, as in this case, we use the notation c7, 


Finally, note the economy in expressing our linear programming problem in terms 
of matrices. You will find that matrix notation gives us an excellent short way to 
express and to simplify problems. 
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Summary of Section 2 


1. The matrix C = AB can be formed if the number of columns of the left-hand 
matrix A equals the number of rows of the right-hand matrix B. 


The product matrix C will have the same number of rows as A, and the same 
number of columns as B, and its elements are given by 


Cy = 1b 15 + Gigbay + °° + GiPny 
a 
= Y ays 
rest 


where n is the number of columns of A. 


2. The commutative law for matrix multiplication is nor true, So, in general, AB 
is not usually equal to BA. However, the multiplicative associative law 


A(BC) = (AB)C 
and the distributive laws 
A(B + C) = AB+ AC 
and 
(B+ C)A=BA+CA 
are true. 
3. The product of a matrix and a column vector is another column vector. This 
enables us to express: 
(i) simultaneous equations in the form 
Ax =b; 
(ii) linear transformations in the form 


r, = Ar,, where r= fa 


4, The transpose of a matrix A is the matrix A? whose columns are the rows of A. 
The main property of transposed matrices is 


(A\Ap..-A,)” = ATAT. |... ATAT, 


End of section exercises 
There are more end of section exercises than usual. This is to ensure that those of 
you who haven't met matrix multiplication before get adequate practice. 


Exercise 15 
Given the matrices 


E22 F 
raat 3 cot ee igs t 
ree alr ens: | } »-[1 § 
=1 2 i ey gly er) 3 
form the following products where possible: 


(i) Ab (ii) AC (iii) ADs (iv) BC (vy) BTC (vi) DEN. 


Exercise 16 


Show that the associative law A(BC) = (AB)C and the two versions of the distributive law, 
A(B + C) = AB + AC and (B + C)A = BA + CA, hold for the matrices 


20nd fey 1 0 = 
A=j-i 2) ii; Be l<t 2 3, Cello 2 1. 
jy ye Ss cba ve a =1 0 2 


(You have already found AB and BA in Exercise 7.) 
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Exercise 17 
Show that the distributive laws stated in the last exercise are true for the matrices 


2 1 1 0 
T 
a-[; 3 + a- [1 '} e-|2 ae 
= 5 1 4 =-1 
Why is the associative law not applicable? 


Exercise 18 
Express the simultaneous equations 
2x, + 3x, + 4x, = -1 
x, — 2x, + 5x5 i} 
x; -2x,= 10 


in the form Ax = b. 


Exercise 19 
Combine the transformations 


sfc l2 Te] [*] [: | 
= id = 
Pl-G aD] = Gl-l a0 
into a single transformation, giving [*| in terms of it 
Ys 


1 
(Solutions to Exercises 15-19 on pp. 57-58) 


3 Change of axes 


3.1 Change of axes in two dimensions (Tape Subsection) 


This subsection is a tape subsection in preparation for the television programme. 
Start the tape now. 
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0 
Given the co-ordinates (x,y) of a point P with respect 


to seme axes, what are the co-ordinates (x',y') of 
this point with respect to some new axes? 


(2) The strategy 


Split. the problem up 


P (uw) 
@ (x',y') 


Put stages | and 2 together = 
Rotation 


(@,y) Translation , (4 v) ————-+ (=', y') 
Stage! Stage 2 
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What are tne co-ordinates of 
P with respect to the 
O'uv axes? 


Answer in the 
back of the cext 


7 


& The rotation matrizx 


= 4 os « 
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use triangle 


O'AQ 


OA + AR 
= OAT as 


+V SIN 


Ps - Rs 


OO ———n— — +a 


+ V sim 


Yoo aw sine + v cos « 


The equivalent matrizxc equation : 


the 
roration 
matrix 
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° Exercise Exercise 2 


P( 2,3) 
@ (x4) 


What are the coordinates of P 
with respect to the O'x'y' axes? 


Answer in the back of text 


Co) Different point of view 


Pp (uv) 


@m ¥ (a) P' has coordinates (u,v), where 


ae = IE 


(») P' has coordinates (uv) where 


uw 


Cos (-«) -Sin (-«) 
sin (-e) cos (-«) 


_ | osn Sind =x 
-sins cose y 


New co ordinates (u,v) of Pare 


Fal 
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‘The Onswer 
in terms of 
x andy 


What are the coordinates of P 
relative to the O'x' y' axes? 


Answer in the 
back of text 
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3.2. The re-orientation problem (Television Subsection) 


Read the following notes before viewing the television programme. 


In the television programme, we solve a three-dimensional problem which is the 
equivalent of the two-dimensional change of axes problem in the tape subsection. TV20 
In Figure 1 below, for instance, we require to find the new coordinates (x', y’,z’) of 

P with respect to an origin O’, given that P has coordinates (x, y,z) with respect to 

the origin O. 


y, Bye) 
(oz) 


) (a,b,c) 


Figure 1 x 


In the television programme, 0 is the planet Earth, P is the planet Saturn, and O' 
is the space probe Voyager. The displacement vector OP of Saturn relative to Earth 
is referred to as s, and the displacement vector OO’ of Voyager relative to Earth is 
referred to as y. Thus from Unit 14 on vector algebra, we know that the 
displacement vector O’P of Saturn relative to Voyager will be s — v. These vectors 
are marked in Figure 2, 


y Paya) 


Figure 2 


At the beginning of the programme, you will see an explicit set of equations which 
determine the new coordinates, written without the use of matrices. You will see 
that these are very long and messy to write down, and have no obvious structure, 
The programme derives a matrix solution to the problem which is (a) more readily 
understandable, (b) easy to derive from first principles, and (c) allows us to see 
some structure in the resulting solution. 


When you watch the programme, you should concentrate on the strategy, rather 
than on the detail of the matrices. These are all given in the post-programme 
notes. 


Now watch the television programme: ‘Applying matrices—the algebra of re- 
orientation’. 


Read the following notes after watching the programme. (They summarize the 
programme, and could help to serve as a substitute for it if you have not been able to 
view it.) 

The problem in the programme is the three-dimensional equivalent of the problem 
on the tape. (See Figure 1 above.) 


P (the planet Saturn) has a known set of coordinates with respect to a set of axes 
Oxyz, centred on Earth. Given another set of axes of known orientation whose 
origin is at a known point O’ (the space probe Voyager), what are the coordinates 
of P with respect to these other axes? 


As before, the problem can be broken down into stages, 
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Stage 1 

Itisa relatively easy matter to get the coordinates of P with respect to axes at O’ 
which are parallel to the original axes, We call these intermediate coordinates 

u, v, W. 


id P (xyz) 
(uvw) 


Figure 3 f 


These can be found by finding the vector OP in terms of the known vectors OP 
and OO’. Using the vector rule of addition, we know that 


OP= OF - 00°. 
In terms of column vectors, this gives us 


x a 
vl=|y]—]b]. (1) 
z ie 


Stage 2 


unr = f ; a P (uv! 
The rotation is not the simple problem that we had in two dimensions, abi 


(xyz) 
The alignment of the O’uvw axes with the O'x'y'z’ axes (see Figure 4) can be done 
in various ways. The way we use is attributed to Euler. The strategy of this 
method is to break the rotation down into three simpler rotations, each of which 
holds one of the coordinate axes fixed. The first rotation we make does not 
actually align any of the axes, but positions the axes in such a way as to ensure 
that the two subsequent rotations will perform the entire re-alignment. The 
advantage of breaking the problem down in this way is that we know a formula 


for two-dimensional rotations. Figure 4 


Imagine that the w-axis in Figure Sis sticking up vertically into the air, and the uv-plane 
is the ground. If a rotation is made, keeping the w-axis fixed, the rotation will not 
change the w coordinate, as the distance of the point above the ground will not change. 
This, and the results from the tape subsection, give us 


u’= ucosa +vsine 
v’ = —usina + v cosa 
w= w, 


In matrix terms, we get the coordinates of P with respect to the axes O'u'v'w’ as 


w cosa sing O}fu 
vl = [ss cosa 0 ‘} (2) 
wi! 0 0 14 Lw. 


If we do a subsequent rotation, this time keeping the u’-axis fixed, and rotating 

through an angle f in the v‘w’-plane, then Figure 6 shows that the u’ coordinate 
will not change, and we get yet another set of coordinates for P, with respect to 
another set of axes, O'u"v"w": 


= u 
= v'cosB + w'sinp 
= —v'sin B + w'cos B. 


ul 
v" 
“ri 


Figure 6 
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This can be expressed in matrix terms as 


u" 1 0 0 ut 
Ca id |: cosp sin | | v’ |. (3) 
w" 0 —sinB cos f| | w’ 

If yet a third rotation is done, keeping w’’ fixed, and rotating about the w’’-axis 
through an angle », then the w”’ coordinate will not change. Similarly to the first 
rotation, we obtain the coordinates of P with respect to yet another set of axes 
Ou" Ww": 


ul’ cos y siny 0] fu’ 
ve’) = -siny cosy 0} |v" |. (4) 
wt 0 0 1 lw 


Suitable choice of the angles «, 8 and y (sometimes called the Eulerian angles) will 
ensure that the O'u"'r""w” axes thus obtained are the required O'x'y’z’ axes, and 
thus that the coordinates (u’"’,v’”’, w’’") are the required coordinates (x',y',2’). How 
these angles are chosen is shown in the programme, and illustrated in Figure 7 


below. (To avoid a confusing figure, the v' and v’’ axes are not shown.) 


Figure 7 


The uv-plane and the x'y’-plane will intersect in a straight line. (Figure 7 shows 
two circular portions of these planes, centred at the origin, intersecting in the line 
marked O'u'.) The three rotations are carried out as follows. 

Step 1 

Rotate about the w-axis until the u-axis is in the line of intersection of the uv- and 
x’'y'-planes. The « is the angle of this rotation. (The new u-axis is called the u’- 
axis.) 

Step 2 

Rotate about the u'-axis until the w'-axis coincides with the z'-axis, Then f is the 
angle of this rotation. 

Step 3 

Rotate about the w"-axis until the u- and v"-axes coincide with the x’- and y'- 
axes respectively. Then » is the angle of this rotation. 


We are now in a position to find an expression for the coordinates (x’, y’,z') in 
terms of the coordinates (u,v, ). We know that 


ina 
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of 
E-f 


and M,, M, and M, are the rotation matrices on pp. 28-29. 


and 


Combining these transformations in the way we did in the last section, and using 
the associative rule, we obtain 


fae u 
| =M,M, * from (3) 
wt 


Lz’ 


u 
=M,M,M, : from (2) 
wv 


Ww 
= mmm : 6) 


w, 


Thus we can multiply M,, M, and M, together to obtain the required relation 
between x’, y’, 2’ and u, v, w in the form 


where 
M=M,M,M,. 
Finally, we know that 


-E- 8 3 


where (x, y,z) are the coordinates of P relative to the original origin 0. So, 
substituting this into the matrix equation (5) above, we obtain the new 
coordinates for the point P in terms of the old coordinates: 


er x-a 
[| -mmiat|) 6] (6) 
2 z=c 

Post-television work 


In the following exercise we ask you to do the final step yourself, thus returning 
full circle to the beginning of the television programme. 


Exercise 4 

Obtain three explicit equations, giving x',y’ and 2’ in terms of x, y,z, the displacement 
coordinates a,b,c, and the Eulerian angles a, p, ), 

[Solution on p. 59] 


End of television programme notes 
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Summary of Section 3 


1. Given that the coordinates of P are (x,y) with respect to some axes Ox, Oy 

with their origin at O, then the coordinates of P with respect to another pair of 
axes O'x’, O'y’ with their origin at (a,b) and inclined at an anticlockwise angle « 
to the original axes (see Figure 8) are given by 


x'|_| cose sina |[x-—a 
y'} L-sina  cosa||y—b 


2. In three dimensions, if we wish to move the origin to the point (a,b,c) and re- 
orientate the axes (see Figure 9) the new coordinates of P are given by 


x’ x—a 
| = maa - 4 
z' z= ¢ 


where a, 8 and y are the Eulerian angles of the rotation and M,, M,, M, are 3 x 3 
rotation matrices each describing a rotation which holds one of the intermediate 
coordinate axes fixed. 


4 Square matrices and their inverses 


4.1 Square matrices 


In this section, we will restrict ourselves to the algebra of square matrices. The 
reason for doing this is that: 

(i) there are important results which are true only for square matrices, and 
(ii) real-life situations often produce square matrix blocks. 


Square matrices of the same size can be both added and multiplied, and thus 
satisfy all the laws of algebra discussed in the first part of this unit. Also, suitably 
sized zeros (0) and negative elements (—A) are already defined. But before we can 
use matrix algebra in a way similar to real number algebra, we still need to find a 
matrix analogous to the number 1, and another matrix analogous to the number 
a~'. This section discusses these matrices, and the resulting matrix algebra. 


(a) The unit (or identity) matrix 


This is a square matrix whose diagonal elements a,,, a,,...4,, are all ones, and 
all other elements are zeros. (The ones are said to lie on the main diagonal, which 
starts in the top left-hand corner.) Examples of unit matrices are 


1 0 0 
[ ‘femfe tol 
0 O 1 
Matrices such as these are referred to as I. The size of I is not normally specified, 
as this is usually clear from the context. 


Unit matrices have the important property that, for any square matrix A, 


The following example demonstrates this property for one particular matrix. 


Example 1 


2 if 
IfA= E 4 then 


and 


x 
Figure 9 
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(b) The inverse of a square matrix A 


A basic property of our real number system is that for each real number a (except 
for 0), an inverse a~' exists so that 


a-ta=1. 


Without this property, it would be impossible to solve equations. This leads us to 
ask whether an inverse matrix A~! exists so that 


AVA=L 


Example 2 


From the tape subsection, we know that if the x- and y-axes are rotated through 
an angle a, then the new coordinates of a point P will be given by 


[7] re ee po | () 


Now, starting with the co-ordinates (x’, y’), we know we can ‘undo’ the process Figure 1 
above by rotating through a clockwise angle «. So we know that 


x|_| cos(—a) — sin(—a) ][ x’ 
~ | =sin(—a) cos (—«) |] _y’ 
_ [cosa —sina |[ x" 
~ [sina cosa |] y’ | 
Substituting this expression for [*| into (1), we get 
x')_[ cosa — sina |[ cosa —sina |[ x’ 
‘| L-sina cosa |] sina cosa || y’ | 
Thus the matrix product 
cosa sina |[ cosa —sina 
—sina cos « |] sin cos & 


leaves all vectors [*] fixed and must therefore be the identity matrix, I. So, if we 


call the original rotation matrix M, then the matrix which ‘undoes’ the rotation 
must be M~. 


Exercise 1 
Given that 
M= [ cosa sine 


cos —sin 
‘ | and oN 2 . ai" 
—sing cos 


sina cos a 


verify by multiplying the matrices together that MN = NM = I. 


Exercise 2 
Given that 


where ad — be ¥ 0, then verify that 
AB=BA=1. 
[Solutions to Exercises 1 and 2 on p. 59] 


From this last exercise, it can be seen that if ad — be is not zero, the 2 x 2 matrix 


1 qd -b 
te 
oe rrsralen ‘| (2) 


MST204 20.4 33 


such that A~*A =I. 


When the inverse of an n x n matrix exists, 


A“A=AA7' =I (3) 


In 2 x 2 matrices where ad — be is zero, then no’ inverse exists; but whenever a 
matrix A has an inverse A~' then the commutative property (3) is true. This will 
be proved more generally in Subsection 4.4, 


Example 3 
2 A 2 1 if 4 -1 2 4 
U: Equati ), the it is = 
sing Equation (2), the inverse ot{ al is {i | but IF dl and 
C | have no inverses, as, in both cases, ad — be = 0, 
Exercise 3 


Using Equation (2), write down the inverses (if possible) of the following matrices: 


4 (ee rf eden gee oe Omm 
Mu [i | e) [-3 | oe [ a} 
[Solution on p. 59) 


One use of the inverse matrix is in the theoretical simplification of matrix 
expressions, and in the theoretical solution of matrix equations, 


If we are given a real number equation 
ax =b, 

then the theoretical solution of this equation can be written 
x=a™"b, 


Similarly, if we are given a set of simultaneous equations, represented by the 
matrix equation 


Ax =b, (4) 
then we can also find a theoretical solution of this equation in the following way. 


Assuming that A~' exists, we can multiply both sides of Equation (4) on the left 
by A“! to give 


A~'Ax = A7'D, 
So 
Ix=A~'b as A~'A =I, 


x=A~'b as Ix =x. 


It is in this kind of way that we use A~' to obtain theoretical simplifications of 
matrix expressions, or theoretical solutions of matrix equations. In many problems, 
this is as far as we need go. However, the next subsection discusses a method for 
finding the numerical value of A~', when this is necessary. 


4.2. A method for finding A~! 


From Subsection 4.1, we know that the theoretical solution of the simultaneous 
equations 


Ax =b (5) 
is 
x=A™'b. 


So, a method for finding A~' would be to solve the equations (5), as the following 
example demonstrates, 
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Example 4 
Find the inverse of the matrix E | by solving the equations 
x, + 0x, =b, E, 
2x, + 3x, =b, E, 
Solution 
Using Gaussian elimination we get 
x, + 0x, =, E, 
E, -2E, Ox, + 3x, = by — 2b, Ex, 


This gives the solution 
E,,+3 x, = 4b. - 9, 
x, =d,. 
In matrix terms, this means that the set of equations 
[ SIe]-bih 
2 aia” Uh 
which is of the form Ax = b, has the solution 
EJ-La sh 
x2} L-3 4JLh 
which is of the form x = A~'b. 


So, by solving the equations Ax = b we find that the inverse of [: 3 is 


2 3 
[4 3] 
sll 
As you know from Unit 9, the elimination process can be done by operations on 
the rows of the matrix Alb. In fact, the whole operation of solving the equations, 


back substitution included, can be done by matrix row operations very similar to 
those used in Gaussian elimination, as the following example demonstrates. 


Example 5 
Once again, we shall find the inverse of a 2 x 2 matrix, namely 


[| 4} 


This time, each step will be written down in matrix form. 


Working Matrix equivalent 


To solue To solve A xX = 


x + 2x, = by E, 1 2\[ x4 
x + 4x2 =by E, 1 4 || x2 


Note that b can be written as Ib. 


1] 


R, 
R, 


(5) 


The elimination 


X, +2x, =b Ey iP Re 
E,-E, Ox, + 2x, =b, —b, Ey R,-R, [0 2][x,| 


The matrix equation (6) is just a matrix form of equations E, and E,,. However, 
with hindsight, we know we could have obtained both these matrices simply by 
subtracting R, from R,, which is the exact equivalent of subtracting E, from E, in 
the equations. 
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To continue with the process of solving the equations: 


Working Matrix equivalent 


First, find x. 


x, +2x,=b, E, i 2)fx.}_ 
Ey, +2 Ox, + x, =+(b, —b,) Ex, Rat? (0 Aix] [— 


Again, the process of finding x, by dividing E,, by two could have been achieved 
by dividing R,, by two in both matrices. 


io 
[ee | 
| agen J 
SS 
[ aes | 
Lome 
PF? 


Now, because we are limited to row operations, we make a slight departure from 
the normal procedure of back substitution. Instead of substituting the value we've 
found for x, into E, to find x,, we subtract a suitable multiple of E,, from E, to 
eliminate x, from the first equation. This is in fact exactly the same as back 
substitution, but it is arranged in a different way. 


Working Matrix equivalent 


Now eliminate x, from E,. 


E, — 2Ex x 1— (62-6) Exe R, — 2Ra, [; Fal -| z= Te] Ri 
X= 3(bz — by) Ex oO 1), -t 4]b.} Rw» 


Once again, it can be seen that the process of subtracting twice E,, from E, could 
have been done directly by subtracting twice R,, from R, in the two matrices. In 
this example the process is now finished. (In a different example, we might still 
have had to divide E,, (or R,,) by a factor to obtain x,.) This gives us that the 


E 1 24. 21 
inverse of [ 1 | is [ } 


4 


The point of this example is that the entire operation of finding the inverse of a 
matrix can be done by row operations simultaneously performed on A and I. 


Starting with the matrices 
A lL 


perform the same linear row operations simultaneously on these 
matrices, If A is ‘reduced’ to I, the same row operations will 
change I to A~'. 


Thus, the example above could be done entirely in terms of matrix row operations 
as follows. 


All 

ee 1 O7R, 

[ 4 0 | R, 

1 2 1 O)R, 

R,-R, [3 2 ke ‘| Ra 
1 2 1 07R, 

R,, +2 E 1 =f i] Ra, 
Sea 0 | es “al R,, 
O 1 |-$ 4)Ry, 

| 


It is advisable if doing these problems by hand to check that the matrix you have 
found really is the inverse, by calculating A~'A or AA“. 


This method for finding the inverse is quite systematic, and like Gaussian 
elimination, it is suitable for implementation on a computer. 
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Example 6 
Find A~', given that 


1 eae 
A=/|2 5 8} 
1 1 7 


Solution 
Start with A I 
fl a 3 1 0 07 R, 
2 5 8 0 1 0} R, 
LI 1 7 0 0 1] R; 
2 3 1 0 OR, 
125) 2 i ro ltRe 
yt ese i dy ee 
2S t 0 OoyR, 
Rent eur i Yall ae 
Ome Gal) 23 gt erie 
rite v2) es 1) %0) joy; 
at 02 |e Again 
Ry+6 10 O 11 -} § JiR 
R,-3Ry [1 2 0] 4 -b -1)R, 
R,-(0 1 OO} -1 ¢ “i Ri, 
oo 0 1|-$ & $JR, 
R,,-®,f1 0 0 # -4 FR» 
ST ce lari €1 -j Ray 
lo 0 u ah 4 $4 Rae 
Finish with I rah! 


(The way we have labelled the rows above keeps track of how they change. Thus 
the first time a row operation is done to R,, for example, it becomes R,,; the 
second time it becomes R,,,. Once you have acquired some practice at doing row 
operations like this, you will probably find it unnecessary to keep labelling the 
rows.) 
Exercise 4 
Use Gaussian elimination-type operations to find the inverses of 
(i) lk RY (ii) [1 3.0 
1 ‘| Ot 
ae) 50: 1 
[Solution on p. 59] 


4.3 The existence of A~! 


As you saw in Subsection 4.1, not all matrices have inverses. If we take a matrix 
which we know has no inverse, and attempt to use the method in Subsection 4.2, 
we shall find out how it breaks down. 


Example 7 
Suppose we attempt to find the inverse of 


Ed 


First get rid of the elements 
under the main diagonal using 
Gaussian elimination, 


Now do the equivalent of 
back substitution. 
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(for which ad — be # 0). Using the method described in Subsection 4.2, we get 
A] I 


ed et OTR: 

F 4 | 0 | R, 

1 fl 1 OTR; 

Ce: we Oilnd | R,, 


At this point, we have to stop, as there is no multiple of R, [0 0] which we can 
subtract from R, to reduce A to I. In fact, if ever a row of zeros occurs in this 
way, We are not going to be able to complete the process to find the inverse. This 
means that the condition for A~* to exist is precisely the same condition that 


Ax=b 


should have a unique solution—namely that the rows of A must be linearly 
independent. 


You might get the feeling by now that everything seems to go wrong if the rows of 
a matrix are linearly dependent. This wouldn't be far from the truth! Since nearly 
every proof in linear algebra assumes that the rows of A are not linearly 
dependent, some terminology is needed to indicate briefly whether this is so. 


A singular matrix is a square matrix whose rows are linearly 
dependent, Conversely (and not surprisingly), 


a non-singular matrix is a square matrix whose rows are linearly 
independent. 


A matrix has an inverse if and only if it is non-singular. 


As you saw in Subsection 4.2, the problem of finding an inverse is effectively the 

same as that of solving a set of simultaneous equations, Thus, the discussion about 

rounding error and ill-conditioning in the unit on simultaneous equations can be 

applied to the problem of finding an inverse. Partial pivoting is used in an effort to _ Partial pivoting is discussed in 
minimize induced ill-conditioning and the subsequent build-up of error. However, Section 3 of Unit 9. 

we are not going to enter further into these questions in this course. 


4.4 Matrix algebra 


The laws we discussed in the earlier sections of this unit, and knowledge of specific 
matrices such as 0, I and A~!, enable us now to do some matrix algebra (as 
compared with matrix arithmetic), So long as we only use known facts about 
matrices we may manipulate the matrix blocks without any regard to their actual 
contents (unless we need to, of course!) 


This subsection consists entirely of examples of matrix algebra, and some exercises 
for you to practise with. 


Example 8 
Multiplying out the expression (A + B)?, we get 
(A +B)? = (A +B)(A + B) 
= A(A + B) + B(A + B) (using the distributive law) 
=A? + AB + BA + B?. 
It is important to note that this is as far as we can go, as in general AB does not 


equal BA, so AB + BA cannot be written as one term, as is the case for real 
numbers. However, there is no ambiguity in referring to AA as A2, etc. 


Exercise 5 

Expand the expression 
(A + B)(A — B). 

[Solution on p. 60} 
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Example 9 
We took as our definition of the inverse matrix that 
ATA=I1. 


Now we want to show that if A~'A = I, then AA~* = I, which was stated without 
proof in Subsection 4.1. The argument can be set out as follows, 


We assume that there is some matrix Q such that 
AQ=I. (7) 


We now want to show that Q = A~'. To do this, we multiply both sides of (7) on 
the left by A7'. 


A~'(AQ) = A“'I, and so 


(A7*A)Q = A~', 
by the associative law and the properties of I. That is, 
Q=at 


So we have proved that for any square matrix A which has an inverse A~' such 
that A~'A = 1, it is also true that 


AAT =I. 
Example 10 
Show for non-singular matrices A and B, that 
(AB)-! = B-'A~!. 
This means that we have to show that 
(B-'A~')(AB) = 1. 
Solution 
(B~'A~!)(AB) = B~'(A~'A)B (using the associative law) 
=B''IB (since A~'A = 1) 
=B'B 
=I (since B~'B = 1). 


So the inverse of AB must be B~'A~'. 


Exercise 6 
Show for non-singular matrices A, B and C that 


(ABC)"! = C-'B-'A-'. 
[Solution on p. 60] 


It'can be shown by induction that for non-singular matrices A,,A2,...,A,; 


(AjAg...A,)7? = AS TAS... AZTAS. 


Exercise 7 (short!) 

(i) What is the inverse of A~!? 
(ii) What is the transpose of AT? 
[Solution on p. 60) 


The last example in this subsection concerns a type of matrix called a symmetric 
matrix. 


A symmetric matrix A is defined to be one for which 
AT=A. 
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An example of a symmetric matrix is 


[i 


and, as you can see, it is symmetric about its main diagonal. 


Example 11 

Show that any matrix Q has the property that QQ’ is symmetric, (Note that the 
product QQ’ can always be formed, as the number of columns of Q is 
automatically equal to the number of rows of Q7.) 


Solution 
We need to show that (QQ7)’ = QQ’. Now, using the fact that (AB)? = BTA 
(see Subsection 2.5), we get 

(QQ7)" = (Q")"Q” 

= QQ’, since (Q’)" = Q from Exercise 7(ii). 

So the result is proved. 
Exercise 8 
Given that A? = I, show that 

(I+ A)(I— A) =0. 


Warning: A? = 1 does not necessarily mean that A = +I. You can check for instance that 
1 


0 | has this property. 


Exercise 9 
Show that for any non-singular square matrix A, 
(AT)" = (A747, 
Exercise 10 
(i) Show that if A is any square matrix, then (A + I) and (A — I) commute, ic. 
(A+ D(A-D=(A-D(A+D. 


(ii) Use this to show that if (A + I) is a non-singular matrix, then (A + 1)~* and (A — I) 
commute, ic. 


(AFD (A-1 =(A-D(AF D7 
[Solutions to Exercises 8-10 on p. 60} 


Summary of Section 4 
1. Definitions 


(i) The unit (or identity) matrix I is a square matrix (whose size depends on the 
context), with ones down the main diagonal and zeros everywhere else. For 
any square matrix A, 


AI=IA=A. 


(ii) A singular matrix is a square matrix whose rows are linearly dependent. A 
non-singular matrix is a square matrix whose rows are linearly independent. 


(iii) The inverse matrix A~! of a non-singular matrix A has the property that 
A‘A=AA7=1L 
Singular matrices (and non-square matrices) have no inverses. 
(iv) A symmetric matrix A is one for which 
=k 
2. Computation of A~* 


The problem of finding A~' is effectively the same as solving the set of equations 
represented by 


Ax=b. 
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We use Gaussian elimination type operations, with partial pivoting. Starting with 
the matrix pair 
All, 


perform the same linear row operations simultaneously on both these matrices. If 
A is reduced to I, then the same row operations will change I to A7!. 


3. A property of inverses 
(A,A2...A,)7) = AS AS. AZ AT. 


End of section exercises 


Exercise 11 


Use Gaussian elimination-type operations to find the inverses (where possible) of the 
following matrices: 


(i) [3 4 (ii) [1 = (iii) [1 0 1 
9 12 3 =1 1 1 0 
0 1 1 
Exercise 12 
If the matrices A, B and X are non-singular and 
B=X"'AX, 
find A in terms of B and X. 
Exercise 13 
Show that if A and B are square matrices such that A and B commute, then: 
(i) A” and B? commute, 
(ii) A7! and B~! commute, 
(iii) A and B~! commute, 
[Solutions to Exercises 11-13 on pp. 60-61} 


5 Introduction to determinants 


In the next few units, we shall be considering some applications of matrices—to 
vibrations of mechanical systems, for example. In addition to the matrix algebra 
considered so far in this unit, some of these applications involve a number 
associated with a square matrix, called its determinant. 


5.1 Definition of 2 x 2 and 3 x 3 determinants 


(a) 2 x 2 determinants 
We have seen in Subsection 4.1 that the 2 x 2 matrix 


a ob 
ea [: | 
has an inverse if and only if the number ad — be is not equal to zero. This number 


ad — be is called the determinant of the matrix A, often denoted by det A. Instead 
of det A, we sometimes write 


a i 


where the straight lines indicate that we are referring to the number ad — be, not 
the matrix. 


Example 1 
1 2 
IfA= [: | then 
1 2 
det A = [ 1 


=(1x 4)—(2x3)=-2 


MST204 20.5 4 


Exercise 1 
Evaluate det A if: 


A 3 1 Ae 1 2 
e a-[ ‘| ©) a-[; Z 


[Solution on p. 61) 


(b) The 3 x 3 case 
A definition of the determinant of a 3 x 3 matrix can be given as follows. If 


SU ts 
A= |az, 932 a3) 


then 


@z2, a 4 ber) ay Aaa 


(1) 


=a, ~ 43 +43 


432 3, 431 G33 43,32 
= Gy 147433 — Ay 423432 — Gy 133 + Gy 2495343, 
oe 43421432 — Gy 342245). 
This way of evaluating det A is often referred to as expanding by the top row. This 
is simply because we write down the first element in the top row (a,,) and 
multiply this by the determinant of the elements left when we cross out the row 
and the column in which a,, lies. 


4rr---8p--— hs 
4a) bert 923 
34 G32 433 


This is simply the determinant 


G22 43 


43233 
We then move to the next element in the top row (a,,), and multiply this by the 
determinant of the remaining elements when we cross out the row and column in 
which a, , lies. 


@y~-—Frr--- 443 
Cea 42 Q23 
43, yn 


This is the determinant 


ber 933 


431 933 
We continue along the top row in this fashion until we reach the last element. As 
you can see from the definition, the signs connecting these terms alternate. A 
general n x n determinant may be defined in a similar way, but for the purposes of 
this section we will confine the discussion mostly to 2 x 2 and 3 x 3 determinants, 


Example 2 
2 1 j 
IfA=0 2 1}, 
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then using Equation (1), 


2 1 3 
det A = |0 2 ih 
3 i 6 
2 1 0 ai Oo 2 
“Pei d- bb deheB il 
= (2 x 11) — (-3) + 3 x -6) 
=22+3-18 
=7, 
Exercise 2 
1 Z 
IfA= |-1 2 3), 
4 1 5 
evaluate (i) det A, (ii) det A’, 
Exercise 3 
1 2 3 
Find det A if A = |4 i 6}. 
7 8 9 


[Solutions to Exercises 2 and 3 on p.61) 


You will see some of the situations in which determinants arise in the next 
subsection, and you will need to be able to evaluate determinants in the unit on 
eigenvalues, 


5.2 Situations in which determinants arise 

Surprisingly, determinants can be used to simplify the formulas occurring in 
various mathematical situations. 

(a) Vector products 

(i) The cross product 


In Unit 14 on vector algebra, you saw that the cross product a x b was given by 
the lengthy expression 


a x b = (a,b, — a3b2)i + (a3b, — 4b3)j + (a,b, — ab, Jk. 
We can re-express the middle term, to get 
a x b = (a:b; — ayb,)i — (a,b3 — a3b,)j + (a,b, — ay, )k 


ay a; ay k Q) 


b, by b, 
This expression can be remembered even more conveniently as a 3 x 3 
determinant: 


ijk 
axb=l|a, a, a,j. (3) 
by b, bs 


Note that we take the liberty of putting i, j and k in the top row, even though they 
aren’t numbers. This is because the resulting form is so easy to memorize. 


Exercise 4 
Use the determinant form of the cross product to evaluate a x b, where 


a=3i-j+4k and b=i-j+k. 
[Solution on p.61] 
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(ii) The triple scalar product 
The triple scalar product a-(b x ¢) which you met in Unit 14 can also be most 
easily remembered as a 3 x 3 determinant. 


Expressing b x ¢ in the style of (2) above: 


b. b. bs]. 
2 1 p+ 
Ca Cy C3 


bxe= 


Lee 2 


b. b. 
a-(b x c)=a,|? by -a,|"! bs a, by, 2 
C2 3 ey oh ey cy 
4 @ 4y 
=|b, by by (4) 
cs Ga 4% 


Again, this is far the simplest way of writing the triple scalar product. 


(b) Areas and volumes 
(i) The area of a triangle in two dimensions 
The area of a triangle in two dimensions can be expressed as a 3 x 3 determinant. 
From Unit 14 you know that 
AB x AC = (\AB\|AC| sin 0)k 
where the triangle ABC is as in Figure 1, 0 is the angle between AB and AC (both 
in the xy-plane), and k is the Cartesian unit vector perpendicular to the xy-plane. 
We can also express the area of the triangle ABC as 
Area of triangle ABC = 4/AB||AC|sin 0 
= HB x AC}, () 
Now, 
AB x AC = (b — a) x (¢— a) 
= (b x c) — (a x c) — (b x a) + (a x a) 
(b x ¢) — (a x c) + (a x b) 
(by the properties of cross products) 


i j kj ji j kl |i j k 
=|b, b, 0}>]a, a, 0}+}a, a, 0 
cy ce, 0} |e, c, 0} |b, b 60 


= ((bye — ¢yb2) — (@yeq ~ ¢442) + (a,b, — bya2))k. 
The coefficient of k in this determinant can be written as the determinant 


a b, 3 
and so from Equation (5), 

1 1 1 
ay by ty 
az bp 
(We have had to write ‘modulus of’ because of the unfortunate similarity between 
the modulus sign and the above way of writing determinants.) 


Area of triangle ABC = modulus of 4 k (6) 


Exercise 5 
Show that the area of the triangle OAB in Figure 1 is 
a, by 
3 ja, | 


[Solution on p. 61} 


“4 MST204 20.5 


(ii) Volumes in three dimensions’ 
The volume of the squashed box in Figure 2 (normally called a parallelepiped) can 
be expressed as a determinant. Suppose the position vectors OA, OB and OC are 
a =a,i+ a,j + a3k, 
b=b,i + b,j + b3k, 
c=ci+ cj + ck 
respectively, and n is a unit vector perpendicular to the base OBDC of the 
parallelepiped. 
We know that the height h of the parallelepiped can be expressed as 
h= |alcos Figure 2 
= |arnl, (7) 


where a, n and @ are as shown in Figure 2. Also, if @ is the angle between the vectors b 
and c in Figure 2, we know from Unit 14 that 


b x ¢ = ((bjlc|sin #)n, (8) 


which is the area of the base of the parallelepiped, times a unit vector n. Now, the 
volume V of our parallelepiped is 


V = height x base 
= ((alcos #) x (|b||e| sin 4) 
= |an| x (|b||c|sin 0) from (7) 
= [a> ({bl[e| sin On)| 
=|a-(b x c)] from (8). 


Thus, using Equation (4), we derive a determinant expression for the volume: 


BR: ais 
V=modulus of |b, 5b, dsj, (9) 
hat CC aL. | 
Exercise 6 
(i) Find the volume of the parallelepiped generated (as in Figure 2) by the position 
vectors 


a=i+k, b=i+] and c=j+ 3k 
(ii) Show for this particular case that 

as(b x c)=c-(a x b) =b-(¢ x a). 
[Solution on p. 61) 


You can see from Exercise 6 that cyclic permutations of a, b and ¢ gave the same 
numerical result. This is fairly obvious if we think about the way we generated the 
volume of the parallelepiped. The base was quite arbitrary. We could have just as 
easily chosen OCEA as the base, and its area would be the magnitude of ¢ x a. So 
long as we are careful not to alter the right-handed order of the vectors, any of the 
expressions in Exercise 6 (ii) will generate the volume of the figure. 


(c) Solutions of simultaneous equations 
You can check that the equations 
Gy 4X, + Qy2Xq = by 
Gy 4X1 + Ay2X2 = by 
have the solution 
yb, — ay2b, 4 4b, — a,b, 


> 2 
922041 — G21942 21442 — 441822 
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provided that the denominators aren’t zero. Writing the equations in matrix form 
Ax = b, the solution can be expressed in determinant form 

A, A 
where A is det A, A, is the determinant of the matrix formed by replacing the first 
column of A by b, and A, is the determinant formed in a similar way by replacing 
the second column of A by b, 


There is a similar formula for the solution of any system Ax = b of n equations in 
n unknowns, known as Cramer’s rule: 
A A. A 
Sele ee tale Cr (10) 
provided that A # 0. Here, A is det A, and A, is the determinant of the matrix 
formed by replacing the kth column of A by b. 


As you will see in Subsection 5.4, it turns out in practice that Gaussian elimination 
is a much more efficient way of solving simultaneous equations. However, 
Cramer's rule provides an interesting theoretical view. For example, expressing the 
solutions in the form of Equation (10), we can rephrase in terms of determinants 
the conditions which we found in Unit 9 for a unique solution, no solution and an 
infinite number of solutions. 


1. If A #0, the equations have a unique solution. 

2. If 4 =0 and one or more of the A,s is non-zero, then the equations are 
inconsistent, and have no solution. 

3. IfA=0 and also A, = A, = A, = 0, then the equations are linearly 
dependent, and will usually have an infinite number of solutions. 


5.3. Properties of determinants 


Determinants have some interesting and useful properties. In this subsection we 
explore some of these for 3 x 3 determinants, using Equation (4) which tells us 
that 

ay a as 
b, by by 
Cy Cy Cy 


=ar(b x c) 


and using some of the properties of vectors which you have already met in Unit 
14, and Subsection 5.2 of this unit. In particular, we shall be using the cyclic 
property which you met in Exercise 6: 


a-(b x c)=b-(c x a) =c-(a x b). (11) 


Property 1: interchanging rows multiplies the determinant by —1. 


To examine the effect of interchanging the second and third row, we look at 
a-(c x b). Now we know that 


exb=-bxe, 

so 
a-(c x b) = —a-(b x c) 
=-A. (12) 

If we change the first two rows, we get 

b-(a x c) =—b-(c x a) 

=-A (from (11)), 

The following exercise completes the possibilities, 
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Exercise 7 

Show that c+ (b x a) = —A. 

[Solution on p. 62) 

So, interchanging any two rows of a determinant changes the sign of the 


‘determinant. We shall refer briefly to this property as P1, and to the following 
properties as P2, P3, etc, 


Property 2: if two rows of a determinant A are identical, then A = 0. 
If row 2 and row 3 are identical, then 
A=a-(b x b) 
=0 (a property of the cross product). 
If row 1 and row 2 are identical, then 
A=a'(axc) 
=c'(axa) (from (11)) 
=0. 
Exercise 8 


Show that if row 1 and row 3 are identical, then A = 0. 
[Solution on p. 62] 


Property 3: if a row of a determinant A is multiplied by k, then the value of this 
determinant is kA. 


We know that multiplying either of the vectors in a dot or cross product by k 
multiplies the result by k. Thus, 
(ka)+(b x ¢) = k(a-(b x ¢)) 
=kA. 
pg en) x ¢) and a+(b x (ke)) are each equal to a-(k(b x ¢)), which is equal 
to kA, 


This result isn’t very important in its own right, but it helps to prove the next 
property, which is very useful. 


Property 4: if one row of A is a linear combination of the others, then A = 0. 
Suppose row | is a linear combination k,b + k,¢ of row 2 and row 3. Then 
A = (k,b + k,c)+(b xc). 
= (k,b)+(b x c) + (kz) (b x c) 
= k,(b-(b x c)) + k,(e-(b x €)) 
=0 (from P2), 
Similarly,.if row 2 is a linear combination k,a + k,¢ of row 1 and row 3, then 
A =a-((k,a + kc) x e) 
=ar((ka) x c) + a+((k,e) x ¢) 
= k,(a-(a x )) + k,(a+(e x e)) (from P3) 
=0 (from P2), 
Exercise 9 
Show that if row 3 is a linear combination of rows 1 and 2, then A = 0. 
[Solution on p, 62] 


The consequences of P4 are wide and will be discussed in Subsection 5.4. 


Property 5: if a row of A is zero, then A =0. 
This is just a statement of P4, if one row were zero times another one. 
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Property 6: if a linear combination of a set of rows of a determinant is added to 
another row not in the set, then the value of the determinant remains unchanged. 


Suppose we add k, times row 1 plus k, times row 2 to row 3, for instance. Then 
a-(b x (¢ + k,a + kb) = a-(b x c) + a+ (b x (ka + kpb)) 
=A+0 (from P4).. 
The other possibilities follow in the same way. 
The last property we discuss in any detail in this section will be proved 
algebraically without using vectors. 


Property 7: det A = det A’. 


Proof 

Aaa Tie bb b b, ob 
oe 2 3 1 3 1 2 

by by by) =a,)? 7 i a 

cy om cy 2 1 3 1 2 


= 44(b3¢3 — bye) — a3(byc, — bse) + a5(b,c2 — be,) 
= 4,(ba¢3 — C9b3) — by (a,¢3 — ¢,45) + c4(a,b3 — bya) 


=o ie ec re Co b, 
3 a3 3 a3 by 
het 
=la, = det A’. 
ay 


Property 7 has an important consequence. Properties 1 to 6 can be applied to 

any determinant. In particular, these properties can be used to describe the effect 
which ROW operations on A” have on the value of det A? (and hence the value of 
det A, using Property 7). Now row operations on AT are column operations on A. 
The consequence of this is that in Properties 1 to 6, the word ‘row’ could be 
replaced by ‘column’, and the statements would still be true. For example, 
rewriting Property 5: if a column of A is zero, then A = 0, 


The following example demonstrates how to use these properties to simplify the 
evaluation of determinants. 

Example 3 

Evaluate det A, where 


1 2. 3 
A=|4 5 6). 
ee 3: |S 


(i) | Using P6 twice, we can subtract 4 times row 1 from row 2 and then 7 times 
tow | from row 3, to obtain 


1 2 3 
det A =| 0 -3 =6|. 
0 a 
Again using P6, we can subtract twice the new row 2 from the new row 3 to 
obtain 
1 2 3 
0 =3 —6/=0 (from P5). 
0 0 0 


(ii) Another way of evaluating this determinant would have used column 
operations. We can use P6 twice, to subtract twice column 1 from column 2 
and then 3 times column 1 from column 3, to obtain 


1 0 0 
det A =|4 -3 —6]. 
7 -6 -12 
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Now, subtracting twice the new column 2 from the new column 3, we obtain 


1 0 0 
dettA=|4 -3 0|=0 (from P5). 
7 -6 0 


There are plenty of other ways of showing that this determinant is zero, but the 
Gaussian elimination technique, which P6 allows us to do, is as good as any. 


A look back to P4 tells us that if the rows of A are linearly dependent then the 
value of the determinant will be zero. The four statements about a square matrix: 
(i) the rows of A are linearly dependent, 

(ii) A is a singular matrix, 

(iii) det A =0, 

(iv) A7! does not exist, 

all say the same thing in different ways. Equally well, the statements: 

(i) _ the rows of A are linearly independent, 

(ii) A is a non-singular matrix, 

(iii) detA #0, 

(iv) A7! exists, 

all say the same thing. We shall find in the next unit how the determinant form of 
this statement is useful. 


There is one other useful property which, for completeness, is stated in this section, 
although not proved. 


Property 8: given two square matrices A and B of the same size, then 
det (AB) = det A det B. 


Finally, although the proofs in this section only show that the properties are true 
for 3 x 3 determinants, they are in fact true for determinants of any size. 


5.4 A general method for the evaluation of determinants 


Any n x n determinant may be evaluated by an extension of the way we evaluated 
a 3 x 3 determinant in Subsection 5.1. However, in general it is much simpler to 
use a procedure analogous to Gaussian elimination. The idea of the method is to 
reduce the matrix to upper triangular form, since the determinant of an upper 
triangular matrix is very easy to evaluate. This is illustrated in the following 
example. 


Example 4 


To evaluate the determinant of an upper triangular matrix A, where 


a obiee¢ 
A= [° d ‘| 
0 oY. 


P7 tells us (among other things) that expanding by columns is just as valid as 
expanding by rows, so we expand ’by the first column to obtain 


de lb oe | 
detA = al, f -off ta 
= adf, 


It can be shown that this principle generalizes to the case of an n x n matrix in 
upper triangular form. We thus have the following result. 


If any square matrix is in upper triangular form (ie. all elements 
below the main diagonal are zeros), then its determinant is just 
the product of its diagonal elements. 
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Now, we know frem P6 that the row operations performed in Gaussian 
elimination without pivoting do not affect the value of the determinant. 
(Remember, the only operation allowed us in Gaussian elimination without 
pivoting is to add or subtract 2 multiple of one row to or from another.) So, to 
find det A, use the operations of Gaussian elimination to reduce A to a matrix U 
in upper triangular form, 


U= 
Then 

det A = detU = 4,13) ...Uy,.- 
Example 5 


Use this technique to evaluate det A, where 


1 1 1 
A=]1 4 | 
2 3 8 


First we do the Gaussian elimination steps 
R,=R,-R, and  R3,=R;—2R). 


We know from P6 that this will not affect the value of the determinant, so 


1 1 1| R, 
detA=|0 3 R,, (Stage 1(a)) 
0 1 Rs, 
Subtracting $R2, from R;, gives 
1 1 1 
det A =|0 3 2 (Stage 1(b)), 
0 o 


At this stage, we can write down the value of det A: 
detA = 1x 3x 1¢ 
= 16. 


For large determinants the amount of work involved in evaluating det A using 
Gaussian elimination is far less than without it. 


Exercise 10 
If 


P27 
cf a | 
229; 


find det A (i) directly, (ii) using Gaussian elimination, 


Exercise 11 
0 2 3 
Use Gaussian elimination to evaluate} 1 3 6}. 
2 2 4 


[Solutions to Exercises 10 and 11 on p. 62] 
Finally, let us take a brief look at Cramer’s method for solving n simultaneous 
equations in n unknowns, 


A, A, 
ENS A ocr 6 i 


where the A,s are the determinants of the form given in Subsection 5.2. To solve a 
set of equations in this way would involve the equivalent of n + 1 lots of Gaussian 
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elimination, whereas solving the equations directly only involves using Gaussian 
elimination once. So, although Cramer's rule is a neat theoretical solution, in 
practice we would hurriedly discard it in favour of Gaussian elimination. 


Summary of Section 5 


A determinant, denoted by det A, is a number associated with a given square 
matrix A. If A is a 2 x 2 matrix, then 


ayy 
det A = *) = ayidz2 — i221. 


a22 


a2 


If A is a 3 x 3 matrix, then 


My, G12 413 
detA =]a2,; 22 Q23 
43, 32 433 
422 a23 3) G23 ard) 422 
=a, - +43 . 
432 G33 ayy 433 34 432 


Larger determinants have similar definitions. 
The most important properties of determinants are: 
(i) if the rows (or columns) of A are linearly dependent, then det A = 0, 


(ii) linear combinations of rows (or columns) of a determinant can be added to 
any other row (or column) without changing its value. 


Using this last property, an n x n determinant det A can be evaluated by reducing 
A to an upper triangular matrix U, using the operations of Gaussian elimination. 
Then 


det A = detU = yy ta... . Upy: 


End of section exercise 


Exercise 12 

Find det A if: 

jks ee ae 

¢ “L-1 4 
1 2 =! 

ji) A=]-1 1 2 
2 =) 1 
Yoel 93 

" 1 2 4 =1 

OES Ge eh 
oO aver 


[Solution on p. 62) 
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6 End of unit test 


As a lot of the terminology in this unit may be new to you, you will find a 
glossary of new terms used in Sections 1-5, and some important results, 
immediately following the test. 


The test is divided into two sections. Section A should not take you longer than 
10 minutes. If you get any of this section wrong, and don't understand why, go 
back and check the appropriate section in the unit before you go on to Section B. 


Section A 
There is only one correct option to each of the following questions. 


1, Which of the following is a 3 x 2 matrix? 
i 1 3 0.2 12 


3 x 
(a) 1 0 3 (b) | 0.3 13 (c) [ ; : il (d) y 
3 0 =1 04 a 


14 
© (@bedef] 


2. Given that 
8 6 xr 
2 1l=ly s|, 
3.4 tal 
the value of s is 
(a) 241 (b) S#14+4 (©) Br+St (d) 1 (e) Bx +2y +32. 
3. ra-[$ ; 6 which ofthe following can be ade toa? 
(a) [14 (b) [2 3 (onisies: 87) fa) Ps aS (e) [*| 
5 24 eed s 
Sama 
1 -1 1), 
Aba 2 os equa to 


2 

« [2 

(a) fo 2 (b) f1 2 (c) (6 6) (d) [; 0 (e) [3 
Ses 5 4 3 6 9 


8 2 4 
5. 2)4 30 2} is equal to 
2 


8 60 4 4 30 2 & 30: 2 
4 20 0 Deon 16 4 10 
(d) |4 6 2 (c) |4 6 2 
2s 4 30° 2) 
i 5 0 2 10 60 
M 
6 Ifx=[x, x, x3] andy =| y, |, then xy equals 
Ys 
x xt 
(a) Deyn + Xay2 + Xays] (b) [x2 yn (c) | x2y2 
Xs Vs XV, 


() [xyr x22 ays] (e) (ee iS al 
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2 13 4 
: ‘ ote fe 
7, Which one of the following can be multiplied on the left by [ 4 56 ar 
(a) [2 1 } (b) [4 ¥ 2) (c) |4 3 6 8 
74 5 4 1 1 4 2 
3 5 <a 
(d) [4 (e) [2 3 
1 4 6 
8% [2 1 3)f38 1  1]isequalto 
3 0 4)i1 5 0 
4 -2 0 
(a) | 19 0 4 (b) 6 5 4 (c) | 19 1 2 
ee] 1 18 -3 1 25-5 3 
25 -S 3 
(d) | 19 25 (e) 6 1 3 
1 =5 3 0 0 
2 3 8: 2 0 
9. A matrix equivalent of the equations 
3x, +2x.+ xy= 4 
Xa =X3 
X1 + 2x2 =x, 
is 
(ey fe 2) tie, 4 () i, oo ef 4 4 
0 0 I}}x.J=]1 ti) 0 1}=}1 
tt 21 beg 1 (ae Wee 1 
(Cy. (xii: 20 ew 4 @ 3 2 t= 1 1)fx, 
x//0 1 -1/=/0 Cathay) Xa 
xyiif 2 =a 0 it! Xs 
NO} (se see’ 1] fx, 4 
0 1 =3] |xz/ =| 0 
1 2 =1}Ix, 0 


10. If Ax = b, and A is non-singular, then 
(a) x=bA™' (by) bx=A™!— (¢) x=A™'b d) X=IA—(€), X= ADAM? 
[Solutions to Section A on p.62) 


Section B 


3 1 K 
1. X=]1 0 0}, evaluate 

LY} 2 3 
(i) X41 Gi) X?+X (iii) X94 x2. 
2. If A is a square matrix such that A? = 0, show that 

A(I+ A)? =A, 

3. If A, Band C are n x n matrices; 
(i) simplify (A + B) — (A + C), 
(ii) show that (A + B)? — A(A + 3B) — B(A — B) = 2(B—A)B. 


4. 
(i) Find the inverse of the matrix 


t 0 © 
Leis 2 0}, 
yord 


(ii) Given that 


1 3 1 
A=]|0 3 Vh 
1 0 2. 
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find the column vector x with the property that 


Bx = Ax +e. 

5. 

(i) If A is a square matrix such that 
A? -A+1=0, 

show that 
A'=I-A, 


(ii) Show that 
(R“'AR)’ = R7'A5R 
where R is any non-singular matrix of the same size as A. Generalize this result, 


6. A square matrix A is said to be orthogonal if AA? = I. 


. +4 | : 
(i) Show that A = [ pe is orthogonal. 
£ - 


(ii) Write down the inverse of this matrix. 


(iii) Show that if A and B are orthogonal matrices of the same size, then AB is also 
orthogonal. 

x 
: | is an n x 1 column vector whose elements are the n variables 
Xn 

XyyeeeyXqy and A is an nx n matrix. 


7. Suppose x = 


Then the matrix 
Q=x"Ax, 
regarded as a function of the variables x,, 
(i) What is the size of Q? 
(ii) Find Q when 


cf = ty 


(iii) Find a matrix B such that x"Bx is the quadratic form 


X», is called a quadratic form, 


xP + 2xpxz + x3. 
(Note; there are several possible solutions to this part.) 
8. (More difficult) 


x vy 
Suppose x = ; and y -| | are column vectors (m x | and n x 1 respectively), 
Xm Yn, 
whose elements are variables, as in Question 7, and that A is an n x m matrix. Show that: 
(i) yTAx = (y7Ax)’, 
(ii), (y — Ax)"(y — Ax) = y"y — 2y"Ax + x7ATAx. 
9, Evaluate the determinants of the following matrices. 
(i) 4 1 (ii) 1 0 -3 (iii) 1 3 0 
—2 8 —f 1 0 0 -1 2 
2 3 =1 3 1 
= 2 1 


(Solutions to Section B on pp. 62-64) 
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Glossary of new terms used in this unit 


Matrix: a rectangular block of numbers. 

Matrices: plural of matrix. 

Column vector: a matrix with only one column. 

An m x n matrix: a matrix with m rows and n columns. 

Element: an individual entry in a matrix, 

aj; the element in the ith row and jth column of the matrix A. 

Linear transformation: a relationship between two column vectors x, and x2 
that can be written in the form x; = Ax). 

0, a zero matrix: a matrix each of whose elements is zero. 

I, a unit matrix: a square matrix with ones down the main diagonal, and 
zeros everywhere else. 

A”! the inverse of A: a matrix with the property A~'A = AA“! =I. 

AT, the transpose of A: a matrix whose columns are the rows of A. 

A singular matrix: a square matrix whose rows are linearly dependent. 

A non-singular matrix: a square matrix whose rows are linearly independent. 
A symmetric matrix: a matrix A with the property AT = A. 

Upper triangular form: refers to a square matrix with zeros below the main 
diagonal. 

An orthogonal matrix: a square matrix A with the property that AAT 
A quadratic form: a single element matrix of the form x"Ax. 


Some matrix results 


Sore 


(AB)" = BTA’. 
(AyAg... An)” 
(AB)! = B-'A7!, 
(Ajaan Ane? AZO ASS AAT. 


Aga ---AZAT. 
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Appendix 


Solutions to the exercises in Section 1 


1. (i) Any matrix with two rows and two columns will 


do—for e: le, . 
lo—for examples) 


ve ae 
(i) |S 6 7 — 8|isa typical example, 
See 203 


1 
2 

(iii) 3/° This is an example of a column vector (or a 
4 


column matrix), 


(ivy) [1 
matrix), 


2. This element would be referred to as pys. (The row is 


always specified before the column.) 


3. By the definition of equality we get 
@ r=4,u=9,2=-3,w=4, 


(i) ax+by=c and dx + ey =f, 
4. (i) A+B 
2 4 Feo ao 
=|4 3 +] 1 21 0 
i 68 6] Oh wed 
2-1 144 340 
=|4+1 341 2+0 
1+0 S41 6+1 
1s; 3] 
=|5 4 2], 
Lead 2] 
(ii) Not possible—as A and C have different sizes, 
(iii) Not possible—as C and x have different sizes. 
2) 7 4 a]_f3 
i woner= (+3) +E )-L+ 
3+a] 
a fi +b] 


implies the equations 


2+1l=a and 4-3=b. 


6. (i) C must be the same size as A and B, so C will be an 


m x n matrix, 


(ii) jj, the element in the ith row and jth column of C, is 
formed by adding the equivalent elements in A and B. So 


cy = ay + by. 


eh 


In this case the zero matrix will be 2 x 2—i.e. i 


0 


a 
b 


4). This is an example of a row vector (or a row 


] 


0 


: 


0 


2 
9% (i) kA=5]3 
4 


10 5] 
(ii) un ents sae 10 -| 10 


2-3 [3 
2-4 =|2 


1 10 5 
2) = ]15 10 |, 
1 20 5 


20 5 
0 5] 
=(|25 10 
30-5. 
10 5 =2 0 
(iii) kA — mB =| 15 10 | — 2 0 
20 =) 2 =2 
12 5 
=|13 10 
18 7 


10. 


(ii) 


i. 


17 7 
=/18 14]. 
25 10 


(i) 2(A + 3B) —4(A—B)= 2A + 6B-4A +4B 


= 10B — 2A. 
X+3(A+X)=2X-A 
4X+3A =2X-—A 
2X = —4A 
X= —-2A. 


(i) (a) a+p-|' 


6 -1 
(b) a-c-[) | 


(c) Not possible. 


wx+ae[™]}+[f]-[2 24] 


ii) 2(P + 2A) + B= 3C; 
2P + 4A + B= 3C; 
P =43C — 4A —B) 


2 S-L “)-E i} 
9) 
= eae 


Solutions to the exercises in Section 2 


1, Using Tables 1 and 2: 
(i) John’s requirements in shop B cost 

(6 x 17) + (3 x 32) + (2 x 33) = 264p, 
(ii) Jane’s requirements in shop A cost 

(10 x 18) + (2 « 30) + (2 x 35) = 310p, 
(iii) Jane's requirements in shop B cost 

(10 x 17) + (2 32) + (2 x 33) = 300p. 
(iv) Joyce's requirements in shop A cost 

(7 x 18) + (2 x 30) + (1 x 35) = 22Ip, 
(v)  Jim’s requirements in shop A cost 

(1S x 18) + (4 x 30) + (3 x 35) = 495p, 
(vi) Jims requirements in shop B cost 

(15 x 17) + (4 x 32) + (3 x 33) = 482p. 
The completed table is in the main text, following Exercise 1.- 


Ce aly eo 
li | 
[4 «3+ (7 = -1) 

~ L(x 3) + (2 x =1) 


Ee 26 
ulin es) 
3. AB has the same number of rows as A, namely 3, and 


the same number of columns as B, namely 2, so AB is a 
3 x 2 matrix. 


(4x 1) + (7 2) 
(1x 1) + (2 x 2) 


4. (i) B will have to have 4 rows to form AB (no 
restrictions on how many columns), so B will be a 4 x n 
matrix. 


(ii) _B will have to have 2 columns to form BA (no 
restrictions on the number of rows), so B will be an m x 2 
matrix. 


5. (i) AB, AC, BC, CA. (BA and CB can't be formed.) 


(ii) ABis a1 x 2 matrix, AC is a 1 x 1 matrix, BC is a 
2 x 1 matrix and CA is a 2 x 2 matrix. 


6. Possible matrix products are BC, BA, CB 
(AC, CA and AB are not possible’) 


BC= [1 -2{ _t] = [5]. 
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Note that AB and BA are not the same. 


8. It is sufficient in these proofs to show that the ijth 
element in the left-hand matrix is the same as the equivalent 
one in the right-hand matrix. We will call such an clement a 
‘typical’ element, 


(i)_A typical element of A(B + C) is obtained by taking the 
ith row of A, 


[a in ve in 
and combining it with the jth column of B + C, 
byt ey 
bay + ey 
: , to obtain 
bay + Gy 


iC + 44) + aialba) + 24) + °° + Aiglbyy + Cy) 


= ¥ aulbyy + cy) 
fen 


" * 
=| ub} + | x auc). 
Tes) wet 
As these expressions are the typical elements of AB and AC, 


A(B + C) = AB + AC. 


(ii) A typical element of (B + C)A is formed by taking the 
ith row of B+ C, 

(bi + eu bin + Cn) 
and combining it with the jth column of A, 


bia + Cia 


ay; 
aay 
.to obtain 
Any, 
(bis + Cia ay + (bin + Cra)ay + °° + (Bin + Cin ny 


= L (bu + cw)ary 


P 

Brea ares) 

As these are typical elements of BA and CA, 
(B+C)A=BA+CA. 


(iii) This is the worst! 


To find a typical element of A(BC), we combine the ith row 
of A, 


[ay js. in 
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with the jth column of BC, 


X bac 
Pest 
giving us the elements 
iy Y binery + Giz Y Dances + °° + din SS bnaces 
fea a1 ket 
au by acu) 
1 Wet 


z 
i o igbgncry 


“ 


=X Y aibertes 
setaet 


== (& Gigban 


Cay 


And this, believe it or not, is a typical element of (AB)C. 


9, 
2 1f2]_fs 
_ [: ‘]-[3} 
10. 
Bot Soe, 2 
at2xat x= 1 
4x, 4 x =-1 
ll. 
2-3 px 1 
7 1 -1{)x:| = a 
3 silo > 
A x= b 
sd 4 = Cry 
= C(Br,) = (CByr, 
= CBiAr,) = (CBA)r, 
=Qr, where Q = CBA. 
=1 ie if 0 
onl Toda bieal Al 
~ -3 -Ss]fi 0 
| 1jl2 -1 
— -13 5 
7. 5 =1) 
13. 
fe Ot 
wf Eo] 
7 1 3 qi) 
14. 


so for these matrices 

(AB)? = BTA’, 
(It is worth noting that AB? would be a 3 x 3 matrix, so 
there would be no hope that (AB)" equalled A’B’,) 


: 2 10/3 7 
wo a-[? IP }-(7] 
Py tape ee 
a ac-[_} ll on :| 
—2 1 4 
-[} 2 | 
(iii) AD can’t be formed. 
(iv) bC can't be formed. 


() b’C=23 ci ; Jee TEST 


1a) aria fo 2 18 
(vi) ver [4 S74 [ - [3 1} 
28 Aa 1. -1 65 
16. From Exercise 7, 
5 8 -1] -—3 Oo 5 
AB=/ 0 6 8 mina | 5 16 -i| 
-1 6 1 5 2 5 
Also, 
2eeeOl eIqin ty ly aid 
aC=|-1 2 1/|0 2 
L3 4 -2J-1 o 2] 
Gat 50: 0] 
m=z) 4 ai, 
L5 8 -7J 
See lee eh 
CA=| 0 2 O]/-1 2 1 
L=t 0 Ls 4 -2J 
f-1 -4 3 
=/-2 4 2], 
4 8 -5) 
po Sie uy i 
BC=/-1 2 3//0 2 0 
et 4 iret oe 2 
f2 4 -3] 
=|-4 4 7 
L2 8 =-1 
Thus: 
0) 0 -1 
(AB)C 2 | 
Qo 2 
4 -3 
A(BC) 4 | 
8 -1 


so for the given matrices, 


(AB)C = A(BC). 
(ii) 
fa io). af 2 92 =2 
A(B+C)=|-1 2 il[- 4 | 
L 4 -2Jb2 4 3 
TG 18; 
=|-2 10 u] 
[-2 14 0. 
fs 8 -1 1 0 0 
AB+AC=| 0 6 i] [2 4 | 
l-7 6 7 5 8 7 
f6 8 =I 
=|-2 10 | 
=? 4 0 


so for these matrices, 
A(B + C) = AB+AC. 


(iii) 
2 22 0 4 
(B+C)A= |-1 4 3) /-1 2 1 
L 4 3aJL3 4 -2 
[-4 -4 8 
=|3 2 -3), 
L9 20 0 
f-3 0 $ -4 
BA+CA=| 5 16 —S 2 
ose ites 8 -5 
f-4 -4 8] 
=| 3 2 -3}, 
L9 2 oJ 


so for these matrices, 
(B+ C)A = BA +CA. 


so 


as 18 3602 
ap+ace[ Br bie | 


= 75 20 
~ L=15 iy 


Now, 


Bi 1 
2 3 FilS 75 (20 
weo? 3 Te J-L2 4 
J 25 0. =15 1 


so the left-handed distributive law is true for these matrices. 
(ii) 


2 1 & 8 aD 
ms [o AE 3 ‘}-[3 0-6 
Sand a Mo 1S 33. 
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so 
re 6 12] 2 3 a 
BA+CA= | 3 Oi OT ily: 6 8 
Ut 15 33) 7 2 30. 
[7 19] 
= |10 6 2 
[18 27 63. 
Now, 
3 1 
wron= i dE | 
19 0. 
[7 a! 19] 
= |10 6 2), 
lis 27 63] 


so the right-handed distributive law is true for these matrices. 


(iii) The associative law does not apply, as BC can't be 
formed (the number of columns of B is not equal to the 
number of rows of C). 


Solutions to the exercises in Section 3 
CEI E)-Epes- 
Pl-Le Sal] 


[2% ably) 


so 


3. 


ele 
[= =IC]-E 


sin 45° |[ 3-2 
cos 45° || 4 — (—5) 
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siny oO) fi 0 0 cosa 
cosy 0) |0 cosB — sinB} | sing 


.66. 
4. cos) sing 0 
M,M,M, = ak cosa 0 


0 0 —sinB — cosf, 0 i) 1 

cosy  sinycosB —_sinysin ff cosx sing =O) 
= |-siny cosycosB cosysinB| |—sine cose 0 
—sinp cos 0 0 1 

cos) cos % — sin) cos B sin x cos) sing + siny cos f cos « sinysin B 

= |-siny cosa — cosy cos f sina —siny sina +cosycosBcosa  cosysinf 

sin B sina —sin B cosa, cos f 

(Alternatively, M,M)M, can be evaluated as M,(M,M,).) y’ =(—siny cosa — cos) cos B sin x)(x — a) 


From Equation (6) on p. 30 we now obtain the equations + (—siny sina + cos cos B cos «)(y — b) 


+ (cosy sin f)(z — c), 
2' = (sin f'sina)(x — a) — (sin Bcos.a)(y — b) 
+ (cos f(z —c). 


= (cosy cos a — siny cos B sina)(x — a) 
+ (cosy sin + siny cos B cos a)(y — b) 
+ (siny sin f)(z — e), 


Solutions to the exercises in Section 4 


(ii) wa-[_ a then a= 1,b = —1,c = —2, and 
i r P d= 2. Thus, 
MN =| 8% sina Jfeosa sina 
“L-=sina — cosa]| sina cos 2 ad — be = (1 x 2)~(—1 x —2)= 
x cos*a + sin?a —cosasina + sinacos«] so the inverse can’t be formed. 
~ | =sinacos.« + cos asin sin? a + cos? 0 1 
[a a (iti) wa-[? jf then a=0, b= t.e= 1d =O, and 
= =I 
OP unles 
a fio) eetlfox a 
Sur el = . 
Ce , F a sls 4 [; | 
cosa —sina cos a sina 
es eae ps ee =| (Note that in this case, A is its own inverse). 
_[  cos?a +sin?« cos asin — sin xcos.a 4) 
~ | sinacos a — cos asin a sin? + cos? AW | 1 
fio 
= lo ‘| =L | 
So MN =NM =I, | 
Dean ep a =F 
bake ry | i 
e) ok ad — be ab + ba 
~ad—beled—de — —cb + da 
yO 
-[5 i]-* 
ad 1 d b b 
ai a 
ne -z{ A ‘| 
_ i . fida—ibe db — bd ; 7 ee 
~ ad—bel-cat+ac ch +ad SoAé -|3 At 
i 0 ~ 
-[5 i] Check either that : ml - ist 
So AB=BA =1. 


F ie a 
2 @ wal? A 


a=2,b=7,c=1,d=6. This gives 


repeal Gees 
~ 5L-1 2y 


I use the formula for A~' on p. 32 with 


Bae aie 


You don't have to do both; remember that A and A~! 
commute. 


(ii) 
A | I 
fl 3 1 5G) 20); ey 
On ey es 0) ak J0l) ees 
Boo 1 ® ior a) oy 
[Cass EO tO) Ol Ry 
O° t =e o 1 Of R 
R,-3R, [0 -9 1 =3) (0) itl! Ry 
i) Rayo t 0 O7 R, 
Dieta 6 1 O| R; 
0 0 9 a 
ji 3 6 1 0 “Oy 
(Meer je i eh|| 
Ry, +28 [0 0 1 wt of ok) Rs 
oe) DO) Ory. 
R,-3R, |0 1 0 & ow a| Ra 
lo oOo 1 zt ote ok] Rs 
i-3R,f1 0 0 * ot of] Ry 
oS 1 0 te ote ome] Ro 
ON foe at wt of oad Rx 
I At 


For safety, you should again check that the A~' you have 
found is correct by finding AA~' or A~'A, and showing it 
comes to I, 
5. 

(A + B)(A — B) = A? + BA — AB - B?, 
(Stop here, as AB # BA in general.) 


6. 
(C™'B-'A~!)(ABC) 
= C~'B'(A~'A)BC (using the associative law) 
Co 'B- IBC (as AW'A =1) 
=C'IC (as B"'B=I) 
=I (asC-'C=D, 


or 


(ABC)~! = (A(BC))~! = (BC)"'A~' = C "BAT! 


7 (i) A. 


As AA~' =I, by definition of an inverse, A must be the 
inverse of A~'. 

(ii) A. 

Let AT = B and (A‘)" = C, and look at a typical element aj 


of A. This becomes the element a, of B (as rows become 
columns), This in turn becomes the element a, of C. 


Thus (A")? =A, 


8. 
(1+ A)U— A) =P + AL—IA— A? 
=P+A-A-I (as A? = 1), 
=0. 


9. We know from Subsection 2.5 that (AB)" = BTA’. 
Letting B= A~', we have 


(AA) = (A“*)PAT, 
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Similarly, 
(Atay = AM(A“4)7 


But (AA~')? and (A~'A)? are both I7 = L. This shows that 
(A~')" is the inverse of A’, which is what we want. 


(A+ D(A 1) =A? +A -AI-P 
=A?-1 


and 


(A-1D(A+1I =A?-IA + AL-P 
=A?-1 
A+D(A—1), 


(ii) From part (i) we have 
(A + 1A —1)=(A— 1A +0). 


Multiply both sides on the left, and also on the right, by 
(A+. 


(A+ D(A +D(A-DA+D 7! 
=(A+ ID (A-1(A+ D(A +1): 


so 


WA = 1A +1)°' = (A +1) 7A DE 
ie, (A-I(AF+ D7 = (A+ Do (A-D. 


3 pale oO] Ry 
R,-3R, [0 0| -3 1} Ray 


A has no inverse. 


(ii) 
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(iii) 
A | I 
pi io. 4 { (0! 407R; 
if et ae 0 1 O}R; 
(Oia 1 © (0 1R 
hea y 1 0 O7R, 
RoR 0 1 —t =I 1 OR, 
1 0 
i) 
1 
1 
1 0 1 1 0 O7R: 
iG 1 1 =1 7 OlR. 
R72 10 O 1 $ -} UR» 
R,-Ry fl 0 0 4 od -$)Ri 
Ry+Rsy ]O 1 0 -$ 4 S/Ry 
o oO 1 4 -t GURy 
1 AS 
Hence 


12, We are given that 
B=X"'AX, 

Multiply both sides on the left by X and on the right by X~!. 
XBX~'=XX"'AXX"'=A, 


Thus, A = XBX~'. 
13. (i) 

A™B! = (BA)! (property of transpose) 
=(AB)' (as A and B commute) 
=B'A' (property of transpose). 

So AT and B’ commute. 
(ii) 

A~'B-' =(BA)~' (property of inverse) 
= (AB)! (as A and B commute) 
=B°'A™! — (property of inverse). 

+ So A~! and B~! commute, 


(iii) 


AB~! =(B-'B)AB™! (multiplying on the left by 1) 


=B'(BA)B"! (using the associative rule) 
=B-'(AB)B-! (as A and B commute) 
=B-'A(BB~') (using the associative rule) 
=B-'A. 


So A and B’! commute. 


Solutions to the exercises in Section 5 


=(3 x 5)—(1 x’ —4) 


=15+4=19. 
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(ii) I J-ux2-exn=0 


2. Using Equation (1); 


1 2 0 
@ |-1 2 3} =i) aay, ; 
2 of 5 ! hues 
= (I x (10 — 3) — (2 x (—5 — 12)) 
=74+3¢=41. 
(ii) 
fi =i! 4 
detaAT=|2 2 1 
3) 4 
ri SR elbs etiam be | 
=1 
E +h +4l i 
= (1% 7) + (1 x 10) + (4 x 6) 
=7+10+24=41. 
3. 
Ha 6 
4 <8) 16 
te eer 
Ss 6] .\4 6 ola 5 
=e ab sll 
=(1 x —3)-(2 x -6)+(3 x —3) 
= -3+12-9=0, 
4. Using Equation (3): 
i> offi» ck 
axb=/3 -1 4 
= 4 
4/3 4). B3 1 
-ft fleet xl 


= 3i+j— 2k. 


5. The triangle OAB is the same as the triangle ABO, so all 
we have to do is to substitute the coordinates of O for those 
of C in Equation (6): 


Area of triangle ABO 


1 1 1 
= {modulus of }a, by | 
a db 0 
- By 2 ON lay» chan) by 
= modulus of {{ ol iB ‘ Co ts i 
=}modutus of | | 
jaa ba} 
6. (i) Using Equation (9): 
1 0 1 
V = modulus of | 1 2 0 
0 1 es 
=(||=7. 
(ii) By part (i), a+(b x ©) = 7. 
0 1 3 1 1 
e*(a x b) 1 0 1j=—- ol* 3 1 
1 2 0 ia 
=(-1 x -1)+ (3x 2)=7, 
2 0 
1 3| Oo 3 
b-(c xa) th a3 =|| |-2 | 
0 1 0 1 1 1 


=(1 x 1)—(2 x -3)=7, 
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Hence 


a+(b xc) =c-(a x b) =bs(c x a). 


c*(a x b) 
=-A 


(property of cross product) 
(using (11)), 


*(a xa) (using (11)) 


(property of cross product). 


9. If row 3 is a linear combination of rows 1 and 2, then 
=k, + kb, and so 
A =a-(b x (kya + kab) 

(b x (kya)) + a+(b x (k2b)) 


= ky(a-(b x a)) + k2(a-(b x b)) (from P3) 
=0 (from P2), 
10. (i) 
1 a 
4 1 3 1 304 
a 4 l= i| | -4 | ry | | 
Bere ie culled rue tue cit ta 
= (1 x 26) — (2 x 22) + (7 x 10) 
= 26 — 44 + 70 = 52. 
(ii) 1 2 2 Es 
3.0064 =2 —20 
=i 2 4 14 
1 2 % 
=|0 -2 =20 
0 0 ~26 
=1x -2x -26=52. 
Ih 
0 2 1 ep 
1 3 Alpes Essential interchange! 
se ote a 2 4| Remember to use P1. 


M 


3 

2 

2 
1 3 6 
-|e 2 3 

4 

3 

2 


o- -8 

1 6 
=-|0 3 

0 0 -2 
=-(1x2x -2)=4 


" 
1 
ie) 


(1 x 3)—(@ x —S)-(1 x =1) 
+10+1 
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i) jl ot 2 a) ft 4 2 3 
ie 32 eel Ge a | 
a os! isi 2 =5 =s 
ool 2 =i) 0) 4 si 
peat 2 3 
o  =4 
=l0 0 =9 3 
0 0 oO 3 

=1x1x 9x3 

= -21, 


Solutions to the end of unit test 
Section A 


1. The answer is (b). Option (c) is wrong because it is a 
2 x 3 matrix. Remember, rows are specified before columns! 


2. The answer is (d). (Equate corresponding coefficients!) 


3. The answer is (d). (Only matrices of the same size can be 
added.) 


4. The answer is (b), (Add corresponding clements,) 
5. The answer is (a). (Multiply each element by 2.) 


6. The answer is (a). (It is the only 1 x 1 matrix among the 
options, for a start!) 


7. The answer is (c), being the only matrix among the 
options with 3 rows, 


8. The answer is (c). 


9. The answer is (c), This can be seen from the fact that the 
equations can be rewritten as 


3x, + 2x, + xy =4 
x, —x; =0 
x, + 2x, — x5 =0, 


10. The answer is (¢), as can be seen by multiplying both 
sides of Ax = b on the left by A~'. 


Section B 
1. (i) 


3 
X+1=]1 
0 


| 


Nee Noe 
wow 
ros 


4 
1 
0 
(ii) 

XP4+X=X(X+1) 


3 1 3] [4 1 
i 0 O;yd 1 
0 2 3) Lo 2 


13 10 21 
2} 3 
Do Gat 
(iii) 


XN? +X? =X(X? + X) 


W 
nou 


W 
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Be 
A(L + A)? = A( + 2A + A?) 
= A(I + 2A) (as A? =0) 
=A+2A? 
ak (as A? =0). 
3. (i) 
(A +B)~(A+C)=A+B-A-C=B-C. 
(ii) 
(A + B)’ — A(A + 3B) — B(A — B) 
= A? + AB + BA + B? — A? — 3AB— BA + B? 
= 2B? —2AB 
= 2(B* — AB) 
= 2B A)B. 
4. (i) 
L | 1 
fd ee ol te 
ay! 0 | 1 | 
woe Oo 0! 8 
fl 0 60 oO © 
nyt 0 | -s 1 i 
Oi Fe Ae Wee 0h 
fio o;1 0 0 
0 4 0 | = 1 | 
0 0 1 | -y+x2 -z 1 
1 | L* 
So the inverse is 
pi o 0 
=x 1 ‘| 
—ytxr  -2z 1 


At this stage we can either use Gaussian elimination (giving 


us x; = 3, x) = 1, xy = 1) or proceed as follows. 
(B— A)x =e; 
x=(B—A)"'s, 
and putting x 1 in the result of part (i), 


Sees 


so I — A is the inverse of A. 
(ii) 
(R~'AR)? = (R~'AR)(R~'AR)(R~AR) 
=R7'A(RR~')A(RR-')AR 
=R'A5R (as RR-! 


(associative rule) 


The general result is 


(R“'ARY=R~-'A"R for any positive integer n. 


6. (i) 
ee a is: 
# -1 \l-4 
J [ ee =e 
Ata et 
(Note, the inverse of an orthogonal matrix is its transpose.) 
(ii) 
A~' =A by part (i) 
J [ ete | 
ee 5 | 
(iii) 
(AB)(AB)’ =(AB)(B7A’) (property of transpose) 
= A(BB")A™ (associative rule) 
=AA™ (because B is an orthogonal matrix) 
=I (because A is an orthogonal matrix), 


‘Thus AB is an orthogonal matrix. 


7. (i) x" is 1 x mand Ax is n x 1. Thus, x7Ax is 1 x 1. 
(In practice, 1 x 1 matrices are usually thought of just as 
numbers.) 


(ii) 


ge bk + dx, 


= [xy(axy + bx2) + x2(ex, + dx2)) 
= [ax} + (b + c)xyx2 + x3). 
(iii) Using part (ii); we must have a= 1,d=1,b +¢=2. 


There are many ways of choosing b and c; for example, 
b=c=1, giving 


8. (i) y"Ax is a matrix with one element, so (y7Ax)" is the 
same one-element matrix, i.e. 


(y"Ax)! = y7Ax, 


(ii) (y — Ax)"(y — Ax) 
=(y? — x7A)(y — Ax) (property of transpose) 
y"y — y"Ax — xTATy + xTATAX 
= y"y — yTAx — (y"Ax)? + xTATAX 
(property of transpose) 
= yy —y™Ax — y"Ax + xATAXx (from part (i)) 
= y"y — 2y"Ax + x7ATAx. 


%» @ he iJ-4 x8 -Gx -2=34 
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(tx -1)-(3 x -5)=14, 


22 


17 


19 


(iii) 


1x -1x -17x -4p2 


= —129, 


